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POST SUNSET BEHAVIOR OF THE 6300 A ATOMIC OXYGEN 

AIRGLOW EMISSION 

I. INTRODUCTION 


Knowledge of the physical, chemical, and dynamical processes of the 
Earth’s upper atmosphere has progressed at a rapid pace during this century. 
The advent of the rocket provided not only the impetus but also, at the same 
time, a means for obtaining some of the observational data required to support 
and verify the theoretical analyses. Knowledge of the "glow" emanating from 
the Earth’s upper atmosphere has shared in this remarkable advancement. 

From the earliest studies it was quite apparent that analyses of the 
twilight glow would provide an opportunity to obtain a wealth of data from 
ground-based observations, because the Sun itself could be used as a differential 
analyzer to establish altitude profiles of the various processes responsible for 
creating the glow. Up to the present, most of the detailed analyses of the 6300 A 
emission in the twilight glow have been concentrated on the pre-dawn twilight 
with little emphasis on the post- sunset phase. This study concentrates on the 
latter situation. 

In thm analysis available observational data were combined with appro- 
priate models and theories in a parametric type of analysis for the purpose of 
establishing a consistent theoretical model that could predict the behavior of 
the intensity of the 6300 A emission during the post-sunset twilight period. 

Analyses of either the pre-dawn or the post-sunset 6300 A twilight glows 
require a well coordinated observational -theoretical approach. Theoretical 
considerations should include: 

1. Chemical reactions, their rate coefficients with their temperature 
dependence, and their efficiencies. 

2. Diffusion and mass transport of the various constituents of the 
atmosphere. 

3. Vertical distributions of both the neutral and charged components of 
the atmosphere and their temporal and spatial variability. 


4 . The thermal structure of both the neutral and charged components. 

5. Solar terrestrial interrelationships, e.g. , sunlight conditions at the 
conjugate point compared with the sunlight conditions at the observational point 
and level of solar and geomagnetic activity at the time of the observation, 

6. Photochemical reactions, their rate coefficients with their tempera- 
ture dependence, and their efficiencies. 

Observational data required either as inputs to the theoretical models of 
the 6300 A emission or for comparison with the calculated results should contain: 

1. Number density-height profiles of both the neutral and charged 
components, 

2. Temperature-height profiles of both the neutral and charged 
components. 

3 . Total integrated emission intensities of the 6300 A emission and the 
nearby background, independently obtained. 

4 . Mass motions, vertically and horizontally, of the neutral and charged 
components. 

5. Photoelectron fluxes and energy spectra entering the atmosphere at 
the observation point. 

Previous theoretical studies and analyses have indicated that dissociative 
photoexcitation of molecular oxygen by solar ultraviolet radiation in the Schumann- 
Runge continuum, dissociative recombination of molecular oxygen and nitric 
oxide ions with thermal electrons, and photoelectron impact excitation of atomic 
oxygen may all produce *D oxygen atoms that contribute to the 6300 A emission 
in the twilight glow. These chemical and photochemical processes are discussed 
in detail in Section n. 

The 6300 A intensity measurements and the equipment used to obtain the 
observational data are discussed in Section HI. 

Constituent number densities at 120 km altitude from Jacchia [1] and 
Reber [2] are combined with analytical temperature profiles after Bates [3] to 
define two of the models of the neutral atmosphere used parametrically in the 
analysis. The third model used is the current Jacchia 1971 model [4] . The 
models and their influence on the computed results are discussed in Sections 
IV and VR. 
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The theoretical model of the 6300 A emission intensity is discussed in 
Section V where the O^D) atoms are assumed to be produced only from disso- 
ciative photoexcitation and dissociative recombination. Consequences of this 
assumption are discussed in Section VIH. 

Solar radiation fluxes in the Schumann-Runge continuum and molecular 
oxygen absorption cross sections from Ackerman [5], Table 14, and solar flux 
values of Friedman and Tousey, as reported by Hinteregger, Hall, and Schmidtke 
[6] with cross sections from Metzger and Cook [7] , Table 13, are used para- 
metrically in this analysis. 

Computational results and comparisons with observational data are 
discussed in detail in Section VII. 

Conclusions are presented in Section VIII with recommendations for future 
research activities outlined in Section IX. 


II. THEORY OF THE 6300 A TWILIGHT GLOW 


A. Sou rce of the Em i ssion 

When the lowest energy level excited state of neutral atomic oxygen 
undergoes a radiative transition to the ground state, a quantum of light is pro- 
duced whose wavelength is proportional to the difference in energy of the levels 
between which the transition takes place. There are three possible radiative 
transitions for the O^D) state: 


, d _6392A_ > _ Spi) 
, d _ 6364A_ 3p_ 


■p 6300A , >p 2 


with respective Einstein coefficients of 1.1 x 10” 6 s -1 , 2.2 x 10“ 3 s -1 , and 6.9 
x 10 “ 3 s" 1 , where the Einstein coefficient is a measure of the probability that a 
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radiative transition will take place between the two levels in 1 s. From the 
foregoing it is apparent that a 6300 A photon is emitted about three times as 
frequently as a 6364 A photon, while the probability of the emission of a 6392 A 
photon is almost negligible. 

Figure 1 shows the low energy portion of the atomic oxygen energy level 
diagram with associated Einstein coefficients, A, and lifetimes, 1/A. Because 
of its relatively long lifetime, the 0( 1 D) state is said to be metastable and, con- 
sequently, it should reach thermal equilibrium with the surrounding neutral con- 
stituents before it undergoes a radiative transition. For this reason the Doppler 
broadening of the integrated 6300 A emission line should provide reliable infor- 
mation about the temperature of the neutral atmosphere in the layer from which 
the greater portion of the radiation is emanating. 

Previous investigators have shown that three mechanisms are capable of 
producing atomic oxygen in the *D excited state. Dissociative photoexcitation of 
molecular oxygen by solar radiation, primarily the S chum ann-R unge continuum, 
can produce both the *D state and the next higher excited state, *3, which in 
turn can produce a *D by a radiative transition in which a 5577 A photon is 
emitted. Dissociative recombinations of molecular oxygen and nitric oxide ions 
with thermal electrons are both postulated to be sources of ^'s with the former 
also capable of producing the *S state; however, the latter reaction, theoreti- 
cally at least, is believed to be incapable of producing the *D state because it 
violates selection rule criteria. Impact excitation of oxygen atoms by photo- 
electrons is also postulated to be a source of the *D and *S states. These three 
production mechanisms are discussed in detail in the remainder of this section. 


B. Dissociative Photoexcitation 

Photodissociation of molecular oxygen by ultraviolet solar photons is the 
primary source of atomic oxygen in the upper atmosphere. The average lifetime 
of the ground state atoms so created is on the order of days. The average 
lifetime of the metastable *D oxygen atoms so produced is about 110 s. The 
Schumann-Runge continuum is the primary energy source for the formation of 
oxygen atoms about 80 km above the Earth’s surface. 

When an ultraviolet photon in the principal Schumann-Runge continuum 
wavelength interval is absorbed by an oxygen molecule, one ground state ( 3 P) 
and one metastable ( 1 D) oxygen atom are produced. The *D then undergoes a 
radiative transition to the ground state with most of the emitted photons having 
a wavelength of 6300 A . 
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Figure 1. Low energy portion of the atomic oxygen energy level diagram [ 8] 
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Brandt [ 9] , on the basis of theoretical calculations, decided that this 
process would account for 25 to 35 kilorayleighs (kR) in the day glow 6300 A 
emission. Noxon [10] , from analysis of ground-based photometer measure- 
ments, concluded that the great variability of the intensity of the 6300 A emission 
appeared to limit the contributions of dissociative photoexcitation to less than a 
few kilorayleighs in the day glow. He set 3 kR as an upper limit on the steady 
contribution from this process. His results also implied that the rate coefficient 
for quenching by O z was nearly 1 x 10 _1 ° cm 3 s” 1 , a result concurred in by 
Dalgamo and Walker [11]. 

Wallace and McElroy [ 12] confirmed previous results which indicated 
that every 0 2 photodissociation by solar flux in the principal Schumann-Runge 

continuum resulted in a *D oxygen atom. 

i \ 

Current observations seqm to place an upper limit of approximately 1 kR 
on the contribution ojf this process to the day glow. 


C. Dissociative Recombination 

Bates [ 13] anA Barbier [ 14] suggested that this process might contribute 
to the observed post-twilight 6300 A emission; Chamberlain [15] showed theo- 
retically that it is probably the primary excitation process in the twilight. 

Brandt [9] proposed that it could account for 15 to 25 kR in the day glow. 
Chamberlain [16] suggested that this process was a major source of the 6300 A 
emission in the night glow. Later, Noxon [10] , in his analysis of ground-based 
photometer measurements in the day and twilight glow, concluded that the great 
variability in the intensity, even under almost identical ionospheric conditions, 
was a strong argument against its being the dominant source. By studying the 
redline intensity variations during the dawn and evening twilights, he showed 
that the observations were quite different from what would be theoretically 
predicted using concurrently measured electron density-height profiles. He 
concluded that the upper limit to contributions from this process would be 
approximately 3 kR in the day glow. 

Despite the consistent picture that is evolving that dissociative recom- 
bination is a major process in the twilight glow, there are still some unresolved 
questions concerning which chemical reactions are involved and which are 
important, the rates at which these reactions proceed, and their efficiencies in 
producing *D oxygen atoms. 
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1. CHEMICAL REACTIONS 

Original concepts involved only the recombination of elect rons with 
molecular oxygen ions, a reaction that can produce the following combinations 
of products: 


o( 3 p) + o( 3 p) 

(1) 

0( 3 P) + 0(*D) 

(2) 

o( 4 d) + o( 4 d) 

(3) 

o( 4 d) + o^s) 

(4) 

o( 3 p) + o( 1 s) 

(5) 


Later, the recombination of electrons with molecular nitric oxide ions process 
was included with the reaction proceeding in either of the two following ways: 

NO + + e — N( 4 S) + O^D) (6) 

or 

NO + + e — N( 2 D) + 0( 3 P) , (7) 

followed by 

N( 2 D) + 0( 3 P) — N( 4 S) + 0( 4 D) . (8) 

Dalgamo and Walker [11] showed that the formation of *D oxygen atoms by 
reaction ( 6) violated the conservation of spin selection rule and also that, from 
an energy standpoint, *S oxygen atoms which could then produce *D oxygen atoms 
through a radiative transition could not be formed. However, Biondi and 
Fiebelman [17] pointed out that there was no direct evidence against the reaction 
to support the theoretical work of Dalgarno and Walker. Wallace and McElroy 
[ 12] supported the work of Dalgarno and Walker [11] by suggesting that quenching 
of N( 2 D) by 0( 3 P) is unlikely and that the most likely agent is 0 2 . 
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Apparently neither reaction (6) or (8) is capable of producing a signifi- 
cant number of *D oxygen atoms; therefore, reaction (6) has not been included 
in this study, and only the molecular oxygen dissociative recombination process 
is considered. 

2. CHEMICAL REACTION RATE COEFFICIENTS 

Rate coefficient values for the above processes are still unresolved. 
Tables 1 and 2 contain values for these rate coefficients that were obtained from 
an extensive literature survey. In this study various combinations of these 
values were used parametrically to establish a consistent set which would pro- 
vide the best agreement between theoretical computations and observations of 
the 6300 A emission intensity in the twilight glow. 

3. PRODUCTION EFFICIENCIES 

Estimates of the efficiency of reactions (1) through (5) for producing *D 
oxygen atoms range from 0. 1 to 1 per each dissociative recombination, being 
almost equally scattered throughout the range of values quoted. Zipf and Fastie 
[46] reported that approximately 10 to 15 percent of the reactions resulted in 
the formation of a *D. Their report was followed chronologically by Tohmatsu, 
Ogawa, and Tsuruta [47] with a value of approximately 16 percent; Peterson 
and Steiger [8] with a value of 10 percent; Brown and Steiger [48] with a value 
of 10 percent; Donahue, Parkinson, Zipf, Doering, Fastie, and Miller [49] with 
a value of 20 percent; Gulledge, Packer, Tilford, and Vanderslice [50] with a 
value of 20 percent; Zipf [51] with a value of 100 percent; Wickwar [36] with a 
value of 50 percent; and Forbes [52] with values of 12 percent at a temperature 
of 865 K, 18 percent at a temperature of 1145 K, and 15±7 percent at a tempera- 
ture of 1000 K. Forbes [52] also concluded that reaction (6) contributed 
approximately 10 percent of 6300 A emission in the night glow. In this study 
efficiencies of both reactions were used parametrically to determine if a "best" 
value could be established. 

4. SOURCE OF MOLECULAR IONS 

Molecular oxygen ions are produced by the charge exchange mechanism 
+ + 

O + O 2 — 1 *■ O + O 2 , (9) 

while nitric oxide ions are produced by the ion-atom interchange mechanism 

0 + + N 2 — NO + + N (10) 
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TABLE 1. VALUES OF BATE COEFFICIENT, a 2 , FOR DISSOCIATIVE 
RECOMBINATION OF O z + PRODUCTION MECHANISM 


0 2 + + e -^*-0* + 0 or 0* + 0* or 0 + 0 


Biondl and Brown [18] 

Sayer l 19] 

Nicolet and Aikin [20] 

Kasner, Rogers, and Biondi [21] 
Danilov [22] 

Whitten and Poppoff [23] 

Ivanov-Kholodny [24] 

Norton, Van Zandt, and Denison [25] 
Bortner [26] 

Swider [27] 

Zipf [28] 

Warke [29] 

DASA [30] 

Kasner and Biondi [ 31 , 32] 

Chan [33] 

G. E. Report [34] 

Biondi [35] 

Wickwar [36] 


2 ( -7) a at 300 K 
4 (-8) at 2500 K 

3 (-8) 

3. 8+1 ( -7) corrected to 1, 7 ( - 7 ) 

3 (-10) T" 1 

4 (-7) -4 (-8) 

1 / 

3 (-7) (300/T) /2 

2 (-7) (300/T) 

2.2 ± 0.5 (-7) (T/300) -0 * 7 
~ 1 (-7) 

1.5 (-7) 

1. 1 (-7) at 300 K 

3 (-7) T" 1,2 

2. 2 ± 0.2 (-7) (T/300) -1 

2.8 ± 0.2 (-7) 

2 ± 2 (-8) (T/300)" 1 to 

4 ± 4 (-8) (T/300)" 1 

1.9 (-7) 

1. 95 (-7) (300/T) 0 * 7 


a. Value in ( ) denotes the order of rate coefficient; e. g. , 2 (-7) means 
2 x 10 -7 . 


9 




TABLE 2. VALUES OF RATE COEFFICIENT, a t , FOR DISSOCIATIVE 
RECOMBINATION OF NO + PRODUCTION MECHANISM 


NO + + e-i^N + 0(*D) 


Nicolet and Aikin [20] 

Gunton and Inn [ 37] 

Lin [38] 

Sy verson, Stern, Shaw, Scheibe 
and Gunton [39] 

Doering and Mahan [40] 

Norton, Van Zandt, and Denison [25] 

Dalgamo and Walker [11] 

Swider [27] 

Saxena [41] 

Weller and Biondi [42, 43] 

Bardsley [44] 

G.E. Report [34] 

Wickwar [36] 


3(-9) a 
1.3 (-6) 

Jr{-9) at 5000 K 

1.3 (-7) at 3000 K 
4 (-7) 

6.7 (-8) (300/T) 

Spin forbidden 

~1 (- 7 ) 

2.8 (-7) 

7.4+ 0.7 (-7) at 200 K 

4.1 (+0.3, -0.2) (-7) at 300 K 

3.1 ± 0.2 (-7) at 400 K 

4.1 (+0.3, -0.2) (-7) (T/300) 

2.6 (-7) at 300 K 

5.0±1 (-7) (T/300) ~ lm 2 ± °* 2 

4.1 (-7) (300/T) 

4.5 (-7) (300/T) 


Meira, Jr. [45] 


a. Value in ( ) denotes the order of rate coefficient; e.g. , 2 (-7) means 
2 x 10 7 . 
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Values of rate coefficients for these reactions are listed in Tables 3 and 4, 
respectively. These values, which were obtained by an extensive literature 
survey, were used parametrically in this study in an effort to establish rr best" 
values. 


Additional sources and sinks of these molecular ions that were obtained 
during the literature survey as well as postulated values of their reaction rate 
coefficients are listed in Table 5 for completeness; however, the production or 
loss mechanisms themselves were determined to be of minor importance in the 
twilight glow, because either the rate coefficient was small in comparison with 
other coefficients or the combination of the reactant number densities and the 
rate coefficient made the mechanism ineffective at the altitudes from which the 
6300 A radiation emanates during the twilight period. 


D. Photoelectron Impact Excitation 

Hanson [53] noticed that photoelectrons emitted above a certain altitude 
could escape from the sunlit atmosphere, travel along the magnetic field lines, 
and be injected into the conjugate atmosphere. Cole [54] then suggested that 
Barbier’s subpolar sheet (pre-dawn enhancement of 6300 A OI emission) could 
be caused by conjugate point photoelectrons causing an increase in the local 
electron temperature, thereby producing an increase in the rate of production 
of 0 1 D atoms by collisions between fast thermal electrons and ground state 
oxygen atoms. This suggestion was contested by Carlson [55] who first reported 
observations of pre-dawn heating of the ionosphere over Arecibo, P.R. , and 
then showed [56] that the enhancement of the 6300 A OI emission had to be due 
to direct impact excitation of atomic oxygen to the *D state by the incident con- 
jugate point photo electrons rather than the fast thermal electrons. This result 
was later supported by the theoretical studies of Fontheim, Beutler, and Nagy 
[86] ; the analysis of simultaneous photometer and Thomson scatter observations 
by Duboin, Lejeune, Petit, and Weill [87]; and the 100 R enhancements reported 
by Cogger and Shepherd [ 88] . If the enhancement of O j D emission had been due 
to heating as proposed by Cole [54], then Duboin, Lejeune, Petit, and Weill [87] 
should have observed only a 3 R enhancement instead of the 40 to 50 R enhance- 
ment they reported. Added evidence in support of the direct impact excitation 
mechanism comes from Wallace and McElroy [12] who concluded from their 
analysis of the measurements from five rocket-borne photometers that this 
mechanism is significant in the day glow, confirming the suggestions of Fournier 
and Nagy [ 89] . It appears as if the major portion of the post-twilight and pre- 
dawn enhancements cannot be attributed to the fast thermal electron mechanism 
proposed by Cole [ 54] . 
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TABLE 3. VALUES OF RATE COEFFICIENT, y 2 , FOR Oa + 0 + 
CHARGE EXCHANGE REACTION 


o 2 + o + _2Vo 2 + + 0 

Bates andNicolet [57,58] 

1("1 2) a 

Dickinson and Sayers [ 59] 

2.5± 0.4 (-11) 

Danilov [22] 

1 (-10) 

Hertzberg [60] 

1 (-9) 

Langstroth and Hasted [61] 

1.8± 0 2 (-12) 

Fite, Rutherford, Snow, and 


Van Lint [62] 

> 1 (-10) 

Norton, Van Zandt, and Denison [25] 

5 (-11) 

Whitten and Poppoff [63] 

2 (-11) 

Swider [27] 

3 (-12) 

Fehsenfeld, Schmeltekopf, and 


Ferguson [64] 

4.0 (-13) 

Donahue [65] 

3.4 (-11) 

Copsey, Smith and Sayers [66] 

2 (-13) 

Wameck [67] 

2 (-13) 

Bohme [68] 

2 (-11) 

Strobel [69] 

1.15 (-11) (600/T) 

Smith and Fouracre [70] 

3. 4 ± 0.5 (-10) T (-°* 48± °* 5 ) 

G.E. Report [34] 

2 (-11) 

Stubbe [71] 

1.16 (-11) (1 + T/1660) 

Bates [72] 

2 (-11) 

Forbes [52] 

2 (-11) increases slightly with 
temperature 

Strobel and McElroy [73] 

1.07 (-9) T -0 * 7 

Meira, Jr. [45] 

2 (-11) 

Biondi [35] 

4 (-11) 


a. Value in ( ) denotes the order of rate coefficient; e.g. , 2 (-7) means 

2x 10“ 7 . 
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TABLE 4. VALUES OF RATE COEFFICIENT, y t , FOR N 2 + 0 + 
ION-ATOM INTERCHANGE REACTION 


o + + n 2 -2Vno + + 

N 

Potter [74] 

~l(-8) a 

Krassovsky [75] 

~1 (- 10 ) 

Bates and Nicolet [ 57, 58] 

~ 1 (-13) 

Hertzberg [60] 

~1 (- 9 ) 

Danilov [22] 

~1 (- 10 ) 

Langstroth and Hasted [61] 

4.7 ± 0.5 (-12) 

Norton, Van Zandt, and Denison [25] 

1 (- 12 ) 

Whitten and Poppoff [63] 

2 (- 12 ) 

Swider [27] 

3 (-13) 

Fehsenfeld, Schmeltekopf, and Ferguson [64] 

3 (-12) 

Copsey, Smith, and Sayers [66] 

2.4 (-12) 

Aquilanti and Volpi [76] 

5.0 ±3 (-12) 

Donahue [65] 

1.8 (- 12 ) 

Waraeck [67] 

2 (- 12 ) 

Paulson, Dale, and Murad [77] 

1.5 (-12) 

Norton [78] 

2 (- 12 ) 

Schmeltekopf, Fehsenfeld, Gilman, and 


Ferguson [79] 

1.8 (- 12 ) 

Smith and Fouracre [70] 

2.4 (-12) 

G. E. Report [34] 

1 (-12 ± 0.5) 

Saxena [41] 

3 (-12) 

Stubbe [71] 

1 (- 12 ) 

Bates [72] 

1.2 (- 12 ) 


a. Value in ( ) denotes the order of rate coefficient; e.g. , 2 (-7) means 
2 x 10 -7 . 
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TABLE 5. OTHER SOURCES AND SINKS AND ASSOCIATED VALUES 
OF THEIR REACTION RATE COEFFICIENTS 


0 2 + + N 2 — NO + 

+ NO 

Fite, Rutherford, Snow, and Van Lint [62] 

Very slow compared with 
charge transfer or ion atom 
interchange 

Ferguson, Fehsenfeld, Goldan, and 
Schmeltekopf [80] 

, <1 (”15) a 

Warneck [67] 

<3 (-15) 

G. E. Report [34] 

1 (-17 ± 2) 

Bates [72] 

<1 (-15) 

0 2 + N + — 0 2 + H 

N 

Danilov [22] 

5.0 (-10) 

Fehsenfeld, Schmeltekopf, and Ferguson 
[64] 

5.0 (-10) 

Copsey, Smith, and Sayers [66] 

4.5 (-10) 

Warneck [67] 

4 (-10) 

Saxena [41] 

5.0 (-10) 

G. E. Report [34] 

5 ± 1.5 (-10) 

O 2 "l* G *“ O 2 

hv 

Nicolet and Swider [81] 

\ 

1-5 (-12) radiative recombi- 
nation can be neglected 

Chan [33] 

2. 8 ± 0.2 (-7) 

G. E. Report [34] 

1 (-12 ± 1) (T/300) -0 * 7 * 0,5 


a. Value in ( ) denotes the order of rate coefficient; e. g. , 2 (-7) means 
2 x 10" 7 . 
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TABLE 5. (Continued) 


NO + + e — NO + hv 

Nicolet and Swider [81] 

1-5 (-12) a radiative recom- 
bination can be neglected 

G. E. Report [34] 

1 ( -12 ± 1) ( T/300) “°* 7 ± °* 5 

O 2 + e -*■ C >2 

+ hv 

Nicolet and Swider [ 81] 

1-5 (-12) radiative recom- 
bination can be neglected 

Chan [33] 

2.8 ± 0.2 (-7) 

G. E. Report [34] 

1 ( -12 ± 1) ( T/300) _0, 7 ± °* 5 

NO + N 2 + — NO + + N 2 

Ferguson, Fehsenfeld, Goldan, and 
Schmeltekopf [80] 

5 (-10) 

Goldan, Schmeltekopf, Fehsenfeld, 
Schiff, and Ferguson [82] 

5 (+1, -3) (-10) 

Wameck [67] 

4.8 (-10) 

G. E. Report [34] 

5.0 ± 1.5 (-10) 

McGowan [83] 

5 (-10) 

Fehsenfeld, Dunkin, and Ferguson [ 84] 

3.4 (-10) ± 30% 

Bates [72] 

5.0 (-10) 

N 2 + + O — * NO + 

+ N 

Norton, Van Zandt, and Denison [25] 

2 (-11) 

Ferguson, Fehsenfeld, Goldan, 
Schmeltekopf, and Schiff [85] 

2.5 ± 1 (-10) 


a. Value in ( ) denotes the order of rate coefficient; e.g. , 2 (-7) means 
2 x 10~ 7 . 
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TABLE 5. (Continued) 


N 2 + + O — NO + + N (continued) 

Wallace and McElroy [12] 

7 (-11)“ 

Saxena [41] 

2.5 (-10) 

G. E. Report [34] 

2.5 ± 0.8 (-10) 

Fehsenfeld, Dunkin, and Ferguson [84] 

1.4 (-10) 

Bates [72] 

2.5 (-10) 

n 2 + + o 2 -* n 2 h 

o 2 + 

Danilov [22] 

1 (-10) 

Fite, Rutherford, Snow, and Van Lint [62] 

2 (-10) 

Fehsenfeld, Schmeltekopf, and 
Ferguson [64] 

1.0 ± 0.5 (-10) 

Wallace and McElroy [12] 

8.6 (-11) 

Aquilanti and Volpi [76] 

6.5 (-11) 

Wameck [67] 

1.1 (-11) 

G. E. Report [34] 

2 ± 1 (-10) 

Bates [72] 

4.7 (-11) 

Meira, Jr. [45] 

i (-io) 

O z + + NO -*• NO + 

+ o 2 

Ferguson, Fehsenfeld, Goldan, and 
Schmeltekopf [80] 

8 (-10) 

Warneck [67] 

7.7 (-10) 


a. Value in ( ) denotes the order of rate coefficient; e. g. , 2 (-7) means 
2 x 10“ 7 . 
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TABLE 5. (Concluded) 


0 2 + + NO — NO + + C 

) 2 (continued) 

Goldan, Schmeltekopf, Fehsenfeld, Schiff, 

8 (+2, -5) (-10) a 

and Ferguson [82] 

G. E. Report [34] 

8.0 ± 2.4 (-10) 

Fehsenfeld, Dunkin, and Ferguson [ 84] 

6.3 (-10) ± 30% 

Bates [72] 

8 (-10) 


a. Value in ( ) denotes the order of rate coefficient; e. g. , 2 (-7) means 
2 x 10" 7 . 


A third source, increased ionospheric recombination due to a lowering 
of the F-region, has been proposed by Noxon and Johanson [90] , Nichol [91] , 
and Pal [ 92] . All but approximately 20 R of the pre-dawn enhancement observed 
over Blue Hill, Massachusetts, could be accounted for by this source when the 
conjugate point was still in darkness, although there were a few enhancements 
that could not be attributed to any of these sources. While most of the 40 R 
enhancement over Hobart, Tasmania, could be attributed to this source, Nichol 
[ 91] also concluded that ( 1) ionization from the exosphere cascading down into 
the F-region after the rapid cooling at sunset, as postulated by Evans [93] , 
was a second possible source and (2) conjugate point photoelectrons were a third 
source, although there was not a consistent correlation between the onset of the 
enhancement and conjugate point sunrise. In view of the recent results of 
Hernandez 1 concerning possible contamination of observations of 6300 A OI 
emission by the 9-3 OH band, which are discussed in detail in Appendix D, the 
possibility that some of these apparently inexplicable enhancements reported by 
Noxon and Johanson [90] as well as the inconsistency between onset time and 
conjugate point sunrise observed by Nichol [91] may be due to increases in the 
intensity of the hydroxyl emission cannot be ruled out, even though the time 
inconsistency does seem to be related to the state of the local ionosphere [91] . 
The amount of enhancement attributed by these investigators to direct impact 
excitation by conjugate point photo electrons is in excellent agreement with the 
results of theoretical studies by Sudworth [ 94] and by Nagy and Banks [95] who 
showed that 13 to 25 R enhancements are possible depending upon the local 
atmospheric conditions. 


1. Hernandez, G.: Contamination of the OI ( 3 P 2 - 1 D 2 ) Emission Line by the 

(9-3) Band of OH x 2 7r in High Resolution Measurements of the Night Sky. 

Unpublished, 1973. 
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The analyses of both Nichol [ 91] and Noxon and Johanson [ 90] suffer 
from a lack of observational data required to positively determine the exact 
amount of the observed intensity due to each of the postulated sources. Neither 
Nichol [91] nor Noxon and Johanson [90] had measurements of the temperatures 
of either the neutral or charged portions of the local atmosphere which are 
required to monitor changing conditions of the environment of the observation 
point nor did they have measurements of the conjugate point photoelectron fluxes. 
The analysis of the pre-dawn enhancement over Arecibo, P. R. , by Wickwar 
[36] was based on the most complete set of observational data since he did have 
the necessaiy electron density data, including data from above the F 2 peak, and 
measurements of conjugate point photoelectron fluxes, although he did not have 
the required temperature measurements either. His analysis showed that when 
the conjugate point solar zenith angle was greater than 105°, the Arecibo night 
glow intensity could be explained by the dissociative recombination process. 
However, the difference between theoretical calculations and observations 
increased with decreasing conjugate point solar zenith angle until a plateau of 
25 to 50 E was reached at approximately at an angle of 95°, supporting Carlson’s 
[96] report that the pre-dawn enhancement over Arecibo began at a conjugate 
point solar zenith angle of approximately 99°. Day-to-day variations in the 
onset time (angle) and intensity were highly correlated with the ionospheric con- 
ditions at the observation point in agreement with the report of Okuda and 
Misawa [97] that when the sunspot number was higher the onset time for pre- 
dawn enhancements over the Takhatta Night Glow Observatory during the period 
from November 1966 to February 1968 was earlier (solar activity effect on 
environment of the observation point) . 

Several investigators have reported on the possibility of geomagnetic 
field control of the location of enhancements. Bennett [98] found a 100 R 
twilight enhancement at solar cycle maximum and a 50 R enhancement at solar 
cycle minimum due to conjugate point photo electrons. He concluded that 
enhancement could not occur at locations with L < 1. 1 because the low minimum 
altitude of the L = 1. 1 shell at 50°W resulted in too many collisions for the 
photoelectrons to be able to propagate from the conjugate point, nor could they 
occur at locations with L > 3. 2 because the electric fields at the plasma-pause 
[ 99] would prevent propagation from the conjugate point. Contrary to this, 

Torr and Torr [ 100] found a 50 R twilight enhancement in winter at the Antarc- 
tica station SANAE, which lies on L = 4.0; Cogger and Shepherd [88] reported 
an enhancement at Saskatoon which lies on L = 4. 2; and Nichol [ 91] found that 
the poleward cutoff of the enhancement was nearer L = 4. 0 than L = 3. 2. A 
solar activity dependent position for the plasmapause could reconcile this 
apparent discrepancy in the poleward cutoff for possible enhancements; however, 
there are still insufficient observational data to verify Bennett's [98] conjecture 
concerning the low latitude cutoff. 
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Huntsville, Alabama, is on approximately the 2. 3 to 2. 4 L-shell, and its 
conjugate point is approximately 47° S and 105° W [101] . Figure 2 shows the 
variation of the conjugate point solar zenith angle as a function of Central Stand- 
ard Time for the night of January 6-7, 1972. If conjugate point photoelectrons 
are present on this night, they should affect the computational results until 
approximately 2233 Central Standard Time, which is conjugate point sundown, 
assuming a screening height of 180 km. 



Figure 2. Conjugate point solar zenith angle as a function of Central 
Standard Time on January 6-7, 1972. 
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E. Quenching 

The metastable *D oxygen atom has a lifetime of approximately 110 s and 
will probably emit a 6300 A photon unless it is collisionally deactivated or 
quenched. The exact process for the quenching, which seems to be a more 
appropriate term for the process, is still undetermined since the *D state is so 
reactive and the emissions are very difficult to observe under laboratory con- 
ditions where rate coefficients can be measured. At the current time, N 2 , 0 2 , 

O, and electrons are postulated to be the primary quenching agents. However, 
the rate coefficient for the 0( 1 D) + O process is relatively slow while the electron 
number densities rarely exceed 10 6 cm -3 at night; therefore, these latter two 
agents do not contribute to the overall quenching process. 

Rate coefficients for the first three of these quenching agents are listed 
in Tables 6, 7, and 8, respectively. Quenching by both N 2 and 0 2 was considered 
to be of importance in the twilight glow; consequently, the reactions listed below 
were included in the computational procedure used in this study with the appro- 
priate rate coefficients used parametrically in order to establish "best" values. 


o(*d) + n 2 — 

»-0( 3 P) + n 2 * 

0( J D) + O z 

— ►o( 3 p) + o 2 * 


( 11 ) 

( 12 ) 


The effects of differences in the quenching rate coefficients on calculated 
volume emission profiles as well as integrated intensities are discussed in 
Section VII. 


III. INSTRUMENTATION AND OBSERVATIONS 


A. Instrumentation 

A general puipose airglow observatory system [123] has been established 
at the Marshall Space Flight Center, Alabama, to study optical emissions of the 
night sky. In this report the system has been used to study post-sunset twilight 
behavior of the 6300 A emission of atomic oxygen. The system consists of an 
airglow observatory, a C-4 vertical incidence ionosonde, and a phase path 
sounder; however, the phase path sounder was not in operation during the period 
of this study. 
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TABLE 6. VALUES OF RATE COEFFICIENT, q*, FOR 

o( 1 d) + n 2 quenching reaction 


0( 1 D) + N 2 -to( 3 P) + N 2 * 

DeMore and Raper [102, 103] 

3(-ll) a 

Hunt [104] 

5 (-15) 

Hunten and McElroy [ 105] 

5 (-11) 

Snelling and Bair [106] 

2.2 ± 1 (-11) 

McGrath and McGarvey [107] 

3 (-11) 

G. E. Report [34] 

5 (-11 ± 0.5) 

Young [108] 

5 (-11) 

Noxon [109] 

5 (-11) 

Peterson and Van Zandt [110] 

5 (-11) - 1 (-10) 

Dalgamo [111] 

Most 0(*D) are quenched by 

n 2 

Forbes [52] 

3 (-11) - 7 (-11) 

5 (-11) - 1 (-10) with 
negligible temperature 
dependence 

Schaeffer, Feldman, and Fastie [112] 

6 (-11) - 7 (-11) 

Donahue, Parkinson, Zipf, Doering, 
Fastie, and Miller [49] 

5 (-11) - 1 (-10) 

Biondi [35] 

4.4 (-11) 

Hernandez [113] 



a. Value in ( ) denotes the order of rate coefficient; e.g, , 2 (-7) means 
2 x 10 -7 . 
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TABLE 7. VALUES OF RATE COEFFICIENT, q 2 , FOR 
O^D) + O z QUENCHING REACTION 


O^D) + 0 2 — ^~0( 3 P) + 0 2 * 

Wallace and Chamberlain [114] 

4 (-12) - 2 (-10) a 

Dalgamo and Walker [11] 

4 (-12) - 2 (-10) 

DeMore andRaper [102,103] 

?(-ll) 

Zipf [115] 

4 (-12) - 2 (10) 

Peterson, Van Zandt, and Norton [ 116] 

> 1 (-U) 

Hunten and McElroy [ 105] 

4 (-15) 

Hunt [104] 

2.5 (-14) 

Warneck and Sullivan [ 117] 

4 (-15) 

Snelling and Bair [ 118] 

2(8) liter mole” 1 s' 1 

G. E. Report [34] 

1 (-12 ± 2) - 1 (-10 ± 2) 

Gilmore, Bauer, and McGowan [119] 

> 1 (-H) 

Zipf [120] 

1 (-H) 

Young, Black, and Slanger [ 108] 

4 (-11) 

Noxon [121] 

6 (-11) 

Biondi [35] 

4.4 (-11) 

Hunt [122] 

5 (-12) 


a. Value in ( ) denotes the order of rate coefficient; e. g. , 2 (-7) means 
2 x 10" 7 . 
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TABLE 8. VALUES OF RATE COEFFICIENT, (fc, FOR 
0( 1 D) + O QUENCHING REACTION 


0( 4 D) + 0( 3 P)-^2 0( 3 P) 

Hunten and McElroy [105] 

1.6 (-13) a 

G. E. Report [34] 

1.3 (-13 ± 0.5) 

Gilmore, Bauer, and McGowan [119] 

< 4 (-13) 


a. Value in ( ) denotes the order of rate coefficient; e.g. , 2 (-7) means 
2 x 10 -7 . 


Three optical measurements are required in this study: (l) the intensity 
of the 6300 A emission, (2) the intensity of the background emission — in this 
study the 6100 A emission is considered to be representative, and (3) the 
Doppler broadening of the 6300 A emission line. Roble, Hays, and Nagy [124] 
have shown that this Doppler broadening of the emission line is due to the 
thermal motions of the metastable 0( 1 D) oxygen atoms and that it is an accurate 
measure of the temperature of the neutral atmosphere at approximately the base 
of the exosphere. The variation of the neutral exospheric temperature during 
the post-sunset twilight period can be monitored by observing the 6300 A Doppler 
temperature of the normal night glow. 

The C-4 vertical incidence ionosonde has been used to obtain electron 
density-height profiles for use in the dissociative recombination portion of the 
theoretical model of the twilight 6300 A emission. Since the C-4 ionosonde is 
only capable of making vertical soundings, only vertical scans of the 6300 A 
emission have been used in this study. 

1. AIRGLOW OBSERVATORY 

The airglow observatory itself consists of an interference filter photom- 
eter after Purdy et al. [125] and a Fabry-Perot interferometer after Roble 
[ 123] . The photometer is capable of measuring the intensity of the night glow 
at four wavelengths — 5300 A , 5577 A , 6100 A , and 6300 A ; however, only the 
latter two measurements have been used in this study. The Fabry-Perot inter- 
ferometer is capable of measuring the Doppler broadening of both the 5577 A 
and 6300 A emissions from the night glow; however, only the latter wavelength 
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emission has been used in this study. Both of these instruments, their design 
considerations, construction, and operating capabilities were discussed in detail 
in Reference 123; however, pertinent portions of that document have been 
included here for completeness. The suitability and capability of the instruments 
for making the required measurements as well as the representativeness of the 
observations have been well documented; therefore, they will not be discussed 
in detail here. Instead, the interested reader is referred to Born and Wolf [126] 
Tolansky [127] ; Jenkins and White [128] ; Turgeon and Shepherd [129] ; Wark 
[130] ; Jarrett, Hoey, and Paffrath [131] ; Jarrett and Hoey [132] ; Hernandez 
[113] ; Biondi and Feibelman [17] ; Roble, Hays, and Nagy [124] ; and Roble 
[123]. 


a. Fabry-Perot Interferometer . This instrument is a conventional, 
pressure scanned, multiple beam interferometer. A simplified block diagram 
of the instrument and its supporting systems is shown in Figure 3. The primary 
components of the instrument are a high quality set of fused quartz etalon plates 
with plane surfaces carefully polished flat to 1/200 of a wavelength and a very 
sensitive photoelectric ally cooled photomultiplier tube. The plates are 15 cm 
in diameter and are separated by a 1 cm spacer which Roble [123] showed to be 
the optimum for measuring the 6300 A emission. Recent observations have 
shown that this may have been an erroneous conclusion. This is discussed in 
detail in Appendix D. The operating parameters of the system are shown in 
Table 9. 

Data from this pulse counting system are displayed on a scaler-timer, 
recorded on a strip chart recorder and an X-Y plotter, and written on a digital 
magnetic tape. 

(l) Data Reduction Procedure. The recorded signal output of the 
6300 A emission line is adequately represented by the following equation which 
is described in detail by Hays and Roble [133] : 


Y(x) 


Io 



+ C 


where 


x is the phase difference; 


I 0 = 5 x 10 5 AS2 R t T P (1 - R) (1 + R)' 1 
u y a T e 7 v 7 


(13) 


(14) 
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Figure 3. Schematic drawing of the Fabry-Perot interferometer and associated 

system components [123]. 










TABLE 9. OPERATING PARAMETERS OF THE FABRY-PEROT 

INTERFEROMETER 


1. E talon Plates 

Effective diameter 
Flatness 

Roughness defect, N^ 

g 

Reflective coatings 

Reflectivity 

Reflective finesse, N„ 

JK 

Spherical defect finesse, N 

D i 

Spacing, t 

2. Objective Lens 

Diameter 

3. Aperture 

Diameter 

Finesse 

4. Instrument Field of View 

5. Interference Filter 

Diameter 
Half- width 
Peak transmission 
Peak wavelength 


5. 25 in. 

X/200 

40 

Five alternate layers of ZnS 
and cryolite 

0.8275 

19.3 

17.3 

1 cm 

48 in. (achromat) 

6 in. 


11/64 in. 
19.6 
0 . 2 ° 


2 in. 
5.2 A 
55% 
6307 A 


TABLE 9. (Concluded) 


6. Photomultiplier (ITT-FW-130) 


Quantum efficiency at 6300 A 

5% 

Photocathode surface 

S-20 

Effective aperture 

1/10 in. 

Dark count (cooled to -15°C) 
(uncooled) 

1 to 2 counts/s 
80 to 90 counts/s 

7. Resolving Power 

350 000 

8. Operating Order 

31 700 

9. Free Spectral Range 

0.198 A (0.5 cm -1 ) 

10. Scanning Gas 

(Dry nitrogen) 

11. Pressure Change for One Order 

1. 76 psi 

12. Overall Instrument Finesse 

10. 00 


where A is the effective collecting area of the interferometer, £2 is the 
instrument field of view, R is the brightness of the source in Rayleighs, t 

y a 

is the transmission coefficient that includes the effects of absorption and 
scattering of the reflective coatings, T T is the transmission of the various 

optical surfaces of the instrument and pre-monochromator, and P g is the 
quantum efficiency of the photomultiplier tube; 


Ao = 7r/2 ; 

A^ = if 1 R n exp[-(n 2 /4) D 2 ] sinc(2nd f /Au) sinc(2nf/Ao) , (15) 
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where R is the reflectivity of the etalon plates, D = 2 7rd /( AcW In 2) with d 

S S 

being the half-width at half-height of the gaussian surface defect function, d^ 

is half-width at half-height of the spherical defect function, A a is the free 
spectral range of the instrument, sine /3= (sin n0) / ( irfi) , and f is half-width at 
half-height of the aperture function; 

n is the number of Fourier transforms required to describe the line 

shape; 


y = 7.274 x l" 12 M _1 o 0 (Atr)" 2 


(16) 


where M is the molecular weight of the emitting atom and <t 0 * s the emission 
wave number; 


T 

n 


is the neutral gas temperature (K) 


y 


and 


C is the continuum background. 


Equations ( 13) and ( 14) were combined and solved for R^ using the following 

values for the various instrument parameters described previously in conjunction 
with the observed intensity of the 6300 A emission above the background 
continuum: 


ASl = 1. 396 x 10" 3 , 


based on an aperture with a diameter of 0. 436 cm, an objective lens with a 
121. 9 cm focal length, a full field of view of 0. 2°, and an effective diameter of 
the etalon plates of 13. 3 cm; 


r = 1 
a 


P =0.05 ; 

e 
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and 



0. 2139 based on the following. 



Optical Component 

Transmission 

(a) 

Glass plate on top of etalon 

0.91 

(b) 

Quartz etalon surface 1 

0.91 

(c) 

Quartz etalon surface 2 

0.91 

(d) 

Glass plate below etalon 

0.91 

(e) 

Objective lens 

0.91 

(f) 

Mirror 

0.91 

(g) 

Collimating lens 

0.91 

(h) 

6300 A filter 

0.55 

(i) 

Condensing lens 

0.91 

0) 

Plexiglass surface 1 

0.91 

(k) 

Plexiglass surface 2 

0.91 


The instrument is designed so that the scanning is accomplished by 
changing the N 2 pressure in the chamber linearly with time, similar to the 
method discussed by Noxon and Johanson [ 133a] . 

b. Photometer . The photometer in the airglow observatory is of the 
turret type designed by Purdy, Megill, and Roach [125] . It has space for six 
interference filters in a wheel that is driven by a variable speed reversible 
stepping motor. This wheel positions an interference filter in front of a photo- 
multiplier tube and then a second wheel positions in succession a dark reference, 
a standard light source, and a sky opening between the interference filter and 
the mirror system that collects the light. Once this procedure is completed, a 
second filter is moved into position and the entire process is repeated. This 
continues uninterrupted during the periods in which measurements are being 
made. 
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The standard light source is a C 14 -activated phosphor [134-136] which 
provides a light of known intensity in the 4000 to 8000 A spectral range. 

E. Marovich of Fritz Peak Observatory calibrated the source against a primary 
standard and its intensity distribution as shown in Table 10 and Figure 4. The 
photometer-telescope system has a 5° field of view and it uses the same scanning 
mirror system as the Fabry-Perot interferometer by having a 45° mirror that 
projects out into, but does not interfere with, the light path from the scanning 
mirror to the Fabry-Perot interferometer. Data from the pulse counting system 
of the photometer are displayed on a scaler-timer, recorded on a strip-chart 
recorder, and written on a digital magnetic tape. 

(1) Data Reduction Procedure. The 6300 A emission line observations 
were combined with the 6100 A continuum observations to determine the intensity 
of the atomic oxygen redline E(6300) in a computer program using the following 
equation developed in Reference 123: 


E(6300) = 


(S -S,)C 
v ss d 7 s 

(S - S ,) 
' sc d 7 


(s - s ,) c 

v cs d 7 c 

(S - S,) 

' cc d 7 


(17) 


where S is the sly signal from the emission line filter, S is the calibra- 
ss J sc 

tion signal from the emission line filter, S is the sky signal from the con- 
es 

tinuum filter, S is the calibration signal from the continuum filter, S , is the 
cc _ _ d 

dark signal of the photomultiplier tube, and C and C are the mean values 

s c 

of the standard light spectral characteristics over the wavelength region of the 
emission line and the continuum as modified by the filters, respectively. 


B. Observational Measurements 


1. INTENSITY 

The Fabry-Perot interferometer and the turret photometer were operated 
during the twilight and post-twilight periods on several days in the winter and 
spring of 1971-1972. The days on which they were operated, the time periods 
covered by the operation, and the solar activity indices existing at the time of 
their operation are listed in Table 11. The photometer system has a 30 A 
bandpass 6300 A filter and measurements are possibly contaminated by other 
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TABLE 10. LOW BRIGHTNESS SOURCE INTENSITY DISTRIBUTION 


Date: March 1, 1972 

Calibration of Phosphor Light: NASA Huntsville Temperature 17°C 


NRD C-3002-a 
Serial D-010 3-71 


Wavelength 

(A) 

Radiance 
(W cm -2 sr -1 A -1 
x 10" 13 ) 

R/A 

Wavelength 

(A) 

Radiance 
(W cm" 2 sr - 1 A _1 
X 10“ 13 ) 

R/A 

3800 

0. 03 

0. 072 

6000 

1.68 

6.38 

3900 

0. 05 

0. 123 

6100 

2.07 

7.99 




6200 

2.34 

9. 18 

4000 

0. 17 

0.430 

6300 

2.49 

9.93 

4100 

0.44 

1. 14 

6400 

2.46 

9.97 

4200 

0. 82 

2. 18 

6500 

2.43 

10.00 

4300 

1. 16 

3. 16 

6600 

2.32 

9.69 

4400 

1.25 

3.48 

6700 

2.02 

8. 57 

4500 

1, 25 

3. 56 

6800 

1.75 

7. 53 

4600 

1. 06 

3.09 

6900 

1. 52 

6.64 

4700 

0. 88 

2. 62 




4800 

0. 73 

2.22 

7000 

1.43 

6.34 

4900 

0. 68 

2. 11 

7100 

1.09 

4.90 




7200 

0.92 

4. 19 

5000 

0. 72 

2. 28 

7300 

0. 73 

3.37 

5100 

0. 68 

2.20 

7400 

0. 53 

2.48 

5200 

0. 68 

2. 24 

7500 

0.41 

1.95 

5300 

0. 67 

2.25 

7600 

0.34 

1.64 

5400 

0. 62 

2. 12 

7700 

0.23 

1. 12 

5500 

0. 61 

2. 12 

7800 

0.24 

1. 18 

5600 

0. 62 

2.20 

7900 

0. 12 

0.600 

5700 

0. 78 

2.81 




5800 

1. 01 

3. 71 

8000 

0.28 

1.42 

5900 

1.32 

4.93 





Notes: 1. The units of radiance are watts per square centimeter per 

steradian per Angstrom. 

2. The precision is to within approximately 10 percent; precision is 
worse toward each end of the wavelength range, especially beyond 
approximately 7500 A . 

3. The accuracy is to within approximately ±10 percent; this is based 
on a partial comparison of the standards. 
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TABLE 11. OPERATION LOG 


Date 

Time Period 

Daily 

10. 7 cm Flux 

81 Day 
Mean Flux 

Kp 

(6.7 h lag) 

Nov. 10, 1971 

1832-0500 

104.8 

109. 0 

0.7 

Nov. 17, 1971 

1720- 1920 

107.9 

106. 1 

0 

Dec. 1, 1971 

1712-1743 

120.3 

110. 5 

1.0 

Jan. 6, 1972 

1726-2400 

106.6 

117. 7 

0.3 

Jan. 7, 1972 

1726-2400 

106.5 

119. 0 

0.4 

Feb. 4, 1972 

1749-1949 

106. 6-105. 6 

119. 1 

2.7-1. 0 

Feb. 7, 1972 

1745-1930 

105.0 

119. 1 

1.0 

Feb. 14, 1972 

1752- 1952 

140.2 

119. 2 

1.0 


emissions; therefore, the intensity values derived using the procedures described 
previously may not be true values for the atomic oxygen redline. A more 
detailed discussion of the contaminating hydroxyl emission is given in Appendix 
D. Each filter in the photometer is calibrated with the low brightness source 
(LBS) immediately after each sky reading and there are only two optical sur- 
faces that must be correctly accounted for in the data reduction scheme, the 
filter and a single mirror, because the scanning mirror system was not used 
for the overhead scans. 


Although the Fabry-Perot interferometer (F-P) has a 5. 2 A half-width 
interference filter that will screen out some more of the contaminating emissions 
from nearby hydroxyl lines, it has 11 optical surfaces that are subject to varia- 
tions with time. Therefore, intensities of the atomic oxygen redline emissions 
were calculated by normalizing the intensities derived from the F-P measure- 
ments using the procedure discussed in Section III. A. 1. a. (l) to the photometer 
intensities during the 2300 to 2330 CST time period where the hydroxyl con- 
tamination should be minimal and then using the normalization factor to correct 
the remaining intensities calculated from the F-P measurements. This is 
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similar to the technique employed by Noxon and Johanson [133a]. Since all of 
the observations used in this analysis were made at the zenith, no corrections 
were made for extinction. It is believed that errors resulting from this pro- 
cedure are rarely more than 10 percent [133a] . Also, it is believed that the 
absolute values of the intensities used in this analysis are accurate to ±10 
percent. 

The observed intensities for the night of January 6, 1972, are shown as 
a function of solar zenith angle in Figure 5. The ±10 percent error bars are 
those due to uncertainties in the calibration of the low brightness source. Errors 
in intensities calculated from the F-P observations were computed using the 
following equation from Reference 133: 


AI 


= (ho±±i_cy/2 


where I 0 and C are defined in Section III. A. 1. a. ( 1) , N is the number of data 
points in each line shape profile, and a' and /3’ are given in Figure 4 of 
Reference 133. In this study the errors ranged from a maximum of 6. 9 R at the 
first observation when the integrated intensity was approximately 250 R to a 
minimum of 1. 2 R at the last two observations when the integrated intensity was 
approximately 18 R. These are all within the ±10 percent errors due to the 
LBS calibration uncertainties. This night was selected for a detailed analysis 
because viewing conditions permitted observations during both the sunset and 
post-sunset twilight periods and the C-4 ionosonde was operational for the 
entire period. Additional observational data are shown in Appendix C. In 
general there is a very rapid decrease in intensity during the first hour after 
sunset followed by a much slower decrease during the second hour and then 
usually a gentle decrease for the remainder of the period. Occasional minor 
perturbations, on the order of a few Rayleighs, are seen in the histories; these 
may be due to either small scale perturbations in the electron densities or to 
wave-like motions that perturb the neutral atmosphere. Determinations of the 
exact causes of these perturbations will require a much more detailed analysis, 
which is not within the scope of this investigation. The absolute values of the 
observed intensities are much more variable from day to day. Once again, a 
more detailed analysis is required to establish the causes of these variations. 
However, in general, it seems that these larger variations are associated with 
the formulation of the Jacchia models of the neutral atmosphere in which both 
the molecular oxygen and nitrogen number densities increase with increasing 
exospheric temperature. This outward expansion of the neutral atmosphere 
during periods of increased solar activity causes an increase in the intensity of 
the 6300 A emission with all other variables, especially the ionosphere, 
remaining constant. 
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Figure 5. 6300 A intensity as a function of solar zenith angle on January 6, 1972, 




2. NEUTRAL TEMPERATURE MEASUREMENTS 

The Doppler broadening of the 6300 A emission line has been shown to be 
a good measure of the temperature of the neutral atmosphere. The procedure 
used to derive temperatures from the Fabry-Perot interferometer measurements 
has been discussed in detail by Hays and Roble [133] . Some typical Fabry-Perot 
interferometer fringes from the twilight period of January 6, 1972, are shown in 
Figures 6 through 13. The figures show intensity counts from the photomultiplier 
tube of the F-P system as a function of pressure in the etalon chamber, with the 
pressure being directly related to wavelength. These are, therefore, scans 
across the 6300 A emission lines, and in theory the half-widths, line widths at 
half-heights, are proportional to the temperature of die emitting particles. In 
Figures 6 and 7, the rapidly decreasing intensity of the combined background 
and line are readily distinguishable. In the remaining figures the decrease in 
intensity appears to be almost negligible. 

One of the assumptions in the data reduction technique discussed in 
References 123 and 133 is that the intensities of both the background and the 
emission line remain constant during the pressure (wavelength) scan. Since 
this is not the case for the first four or five fringes in the twilight period, the 
varying intensities had to be corrected before the data could be reduced. Both 
log and exponential curves were fit through the intensity counts from the F-P 
system during those portions of the pressure scans when only the background 
was being measured. The curve which fit the data best was then used to correct 
for the varying background intensity. Each fringe that showed a variable back- 
ground was reconstructed by subtracting from each intensity data point the 
intensity value computed from the curve fit through the background intensity 
values only. Arbitrarily then, a value of 200 was added to each data point to 
preclude the possibility of obtaining negative intensity values. Figures 14 
through 17 show the results of this reconstruction. These four figures 
correspond to Figures 6 through 9, respectively, i. e. , the immediate post- 
sunset portion of the January 6, 1972, twilight period. In the immediate post- 
sunset period, the 6300 A line is so intense that pressure scans can be accom- 
plished in less than approximately 2 min; therefore, the varying intensity of 
the line itself during the scan period is so small that no additional correction 
is required. 

These "raw" data were then reduced using a computerized version of 
the technique discussed in References 123 and 133. The derived temperatures 
are shown as a function of solar zenith angle in Figure 18. 
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Figure 6. 


Figure 7. 



Fabry-Perot interferometer observations of the 6300 A emission 
line shapes at 1728 and 1732 CST on January 6, 1972. 



Fabry-Perot interferometer observations of the 6300 A emission 
line shapes at 1740 and 1742 CST on January 6, 1972. 
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Figure 8. 



Fabry-Perot interferometer observations of the 6300 A emission 
line shapes at 1751 and 1756 CST on January 6, 1972. 



Figure 9. Fabry-Perot interferometer observations of the 6300 A emission 
line shapes at 1805 and 1811 CST on January 6, 1972. 




Figure 11. Fabry-Perot interferometer observations of the 6300 A emission 
line shapes at 1845 and 1855 CST on January 6, 1972. 
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ETALON PRESSURE 


Figure 12. Fabry-Perot interferometer observations of the 6300 A emission 
line shapes at 1903 and 1910 CST on January 6, 1972. 



ETALON PRESSURE 

Figure 13. Fabry-Perot interferometer observations of the 6300 A emission 
line shapes at 1926 and 1940 CST on January 6, 1972. 
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Figure 14. 6300 A emission line shapes at 1728 and 1732 CST on 
January 6, 1972, after correcting for and removal of the 
variable intensity of the background. 



Figure 15. 6300 A emission line shapes at 1740 and 1742 CST on 
January 6, 1972, after correcting for and removal of the 
variable intensity of the background. 
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Figure 16. 6aioO A emission line shapes at 1751 and 1756 CST on 
January 6, 1972, after correcting for and removal of the 
variable intensity of the background. 



Figure 17. 6300 A emission line shapes at 1805 and 1811 CST on 
January 6, 1972, after correcting for and removal of the 
variable intensity of the background. 
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The error bars shown on the figure were calculated from the following 
equation from Reference 133: 



where Iq and C are defined in Section HI. A. 1. a. ( l) , N is the number of data 
points in the line shape profile, and a and (3 are given in Figure 3 of 
Reference 133. 

3. ELECTRON DENSITY MEASUREMENTS 

The electron density-height profiles used in this theoretical analysis 
were obtained from C-4 ionosonde measurements. True height analysis was 
accomplished using a computer program based on the technique proposed by 
Doupnik and Schmerling [137] . The electron density-tme height profiles above 
the F 2 peak were computed using a constant scale height two times that of the 
atomic oxygen scale height at the F 2 peak with a kinetic temperature of 1500 K. 

Figure 19 is a plot of contours of constant electron density versus solar 
zenith angle for the approximate 7 h period following sunset. Both the height of 
the F 2 peak and the maximum electron density decrease relatively smoothly with 
time, and no large scale motions of the charged particles are apparent in this 
record. 


During the course of this investigation much effort was spent in estab- 
lishing the electron density-true height profiles. Several different reduction 
techniques were compared using known electron density-true height profiles as 
input data. Various starting techniques, amounts of input data, location of 
input data, and types of input data (X-ray only, O-ray only, X- and O-ray data) 
were used in the available data reduction programs to determine which output 
results most closely matched the known profiles. The electron density-true 
height profiles referenced above are believed to be most representative of the 
ionospheric conditions existing over Huntsville, Alabama, on the night of 
January 6, 1972. They were determined by using the combination of starting 
point, amount, type, and location of input data established during the study as 
the one combination that produced output results which most nearly matched the 
known electron density-true height profiles. Although these results are satis- 
factory insofar as this analysis is concerned, this area requires a concerted 
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Figure 19. Contours of constant electron number density. 


research effort. Until such time as this additional research is concluded, 
electron density-true height profiles, other than possibly those reduced by the 
group at NOAA, Boulder, Colorado, should not be used without additional 
verification. 
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IV. NEUTRAL ATMOSPHERE 


A. Models 

Three models of the neutral upper atmosphere of the Earth were used in 
this analysis. The original J acchia 1964 model and the same model with the 
120 km boundary conditions recommended by Reber 2 both used the analytical 
representation of the neutral gas temperature-height profile developed by Bates 
[3] and modified by Walker [138] to calculate constituent number density-height 
profiles. The third model used was the Jacchia 1971 model which was used in 
its original form. Table 12 contains the boundary condition values used as 
inputs to the J acchia 1964 model atmosphere. Corresponding number densities 
at 120 km altitude for the constituents in the Jacchia 1971 model atmosphere are 
variable and dependent upon the exospheric temperature. Differences in the 
computed 6300 A emission integrated intensities and volume emission profiles 
only from the differences inherent in these models are discussed in detail in 
Section VII. 

TABLE 12. NEUTRAL CONSTITUENT NUMBER DENSITIES AT 120 km 


N 2 

o 2 

O 

Source 

3. 7131 x 10 11 cm- 3 

6. 0645 x 10 10 cm -3 

9.4607 x 10 10 cm” 3 

Reber [2] 

4. 0 x 10 11 cm” 3 

7. 5 x 10 10 cm” 3 

7.6 x 10 10 cm" 3 

Jacchia [1] 


B. Columnar Number Densities 

The production of *D oxygen atoms and the subsequent emission of the 
6300 A photons are dependent upon the number densities of the neutral con- 
stituents as well as their total densities in a unit column between the altitude of 
the emitting atoms and the source of the agent causing the excitation, the Sun 
in the case of dissociative photo excitation and the magnetic conjugate point in 
the case of photoelectrons. 


2. Private communication, 1973. 
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To calculate the total columnar number densities for use in the dissocia- 
tive photoexcitation calculations, information is required on the constituent 
number density variation with height in the atmosphere. All major atmospheric 
constituents are believed to be in diffusive equilibrium above approximately 
120 km altitude; therefore, at any altitude, z, above 120 km the constituent 
number density, n.(z) , can be calculated from 


n j( z ) = n i ( z 0 ) 


T(zp) 

T(z) 


z 

exp - f 
z 0 


dz 

Hj(z) 


( 18 ) 


where z 0 is a reference height in the diffusive equilibrium region of the atmos- 
phere, taken as 120 km in this study; T is the kinetic temperature; and H.(z) , 
the scale height of the ith constituent at altitude, z, is defined as 


H i (z) = RT(z)/M.g(z) 


( 19 ) 


where R is the universal gas constant, is the molecular mass of the ith 

constituent, and g( z) , the acceleration due to gravity at the height z in 
kilometers, is defined as 

g( 2 ) = 9-30665xl0-»( 63 6 5 3 6 5 ^ ) . (20) 

From the perfect gas law and the hydrostatic equation, the total number 
density in a unit vertical column above a height, z, in an isothermal atmosphere 
is 


N.(z) = n t (z) H^z) 


( 21 ) 


For nonvertical incidence during twilight observations, the zenith angle 
of the Sun must be known. For geographic latitude, 0 , local time ( angle from 
noon) , 'k , and solar declination, 5 , the zenith angle, x » is obtained from 
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cos x 


cose cos6 cos^ + sine sin6 


( 22 ) 


During the twilight period studied here, the zenith angle ranged from 
97. 5° to 166°. When the zenith angle, x » is 60° or less, the curvature of the 
atmosphere may be neglected and the constituent total columnar number density 
above an altitude, z, is approximated by 


N.(z) * n.(z) H.(z) sec x 


(23) 


However, this is not the case during the twilight period and the curvature 
of the atmosphere must be included in the calculations. Chapman [139, 140] 
first investigated the problem of calculating total columnar number densities in 
a slant path through a spherical isothermal atmosphere. In an isothermal 
atmosphere the scale height, , for a given constituent is constant and indepen- 
dent of height. If the constituent is in hydrostatic equilibrium, the total columnar 
number density above an altitude z , can be calculated from 



where x(l) = (R + z)/H.(z) and R is the Earth’s radius. The Chapman 

Jl( X Hi 

function, generally denoted as Ch(x,x), is the quantity between the braces. It 
was developed on the assumption of an isothermal atmosphere, a situation 
which exists only above about 250 km in the Earth’s atmosphere. Since the 
height profile of the constituent number density, n.(z), is almost exponential, 

the major contribution to the total column density comes from a narrow altitude 
range near z over which the temperature is nearly constant. The total column 
density is relatively insensitive to the temperature-height profile; therefore, 
the constituent scale height at the observational altitude can be used in the 
Chapman function. A number of approximations to the Chapman function have 
been developed; however, the one used in this study was introduced by Rishbeth 
and Garriott [ 141] : 
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Ch(x, x ) = [ir/2 x(l) sin x ] ^ ejq) cos 2 x J |l±erfj^^ cos^J^j 

(25) 

where the symbols are as previously defined and the plus sign is used for zenith 
a ngl es >90° while the minus sign applies to zenith angles <90°. 


C. Exospheric Temperatures 

The exospheric temperatures required to establish Ihe constituent num- 
ber density height profiles in the atmospheric models used in the analysis were 
derived from the analysis of the Doppler broadened emission line profiles 
measured with the Fabry- Perot interferometer after the varying background 
correction had been made. 


V. THEORETICAL APPROACH 


The theoretical model established to study the 6300 A emission intensities 
during the twilight period contains two production sources of *D oxygen atoms. 


A. Photodissociation 

The number of *D oxygen atoms produced at each altitude level by solar 
photons in the principal Schumann-Runge continuum was calculated by 



1775 

Y, <t(X,0 2 ) x FLUX(X) x n(0 2 ) 
\=1125 Z Z 


(26) 


where o(\,0 2 ) is the molecular oxygen absorption cross section for radiation 

of wavelength and FLUX( \) is the photon flux of wavelength \ at height z. 

z 

The photon flux of wavelength X at each altitude z was calculated from 
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» (27) 


FLUX(^) z 


= FLUXCA)^ X EXP 



-N(l) z 


x 



where FLUX ( X) is the flux of solar photons of wavelength X at the top of the 

atmosphere, 500 km in this study; N(I) is the column density of the constituent 

z 

I above level z; and a (A., I) is the absorption cross section of the constituent I 
for photons of wavelength X . The column density of the constituent I is cal- 
culated from 


N(I) = n(I) x H(l) x CH(I) , (28) 

Li Z La Li 

where n(I) is the number density of the constituent I at altitude z; H(I) 
z z 

is its scale height at altitude z computed from 


RT 

H W Z - 5555- 

z 


(29) 


and CH(I) , the Chapman function for constituent I which is used to correct 

for the atmosphere at the large zenith angles during the twilight period, is 
calculated from equation (25) . Each photodissociation was assumed to produce 
one *D oxygen atom. 


B. Dissociative Recombination 

Dissociative recombination of both 0 2 + and NO + with thermal electrons, 
as discussed in Section II, where both of these ions are produced as a result of 
collisions between atomic oxygen ions and oxygen and nitrogen molecules [equa- 
tions (9) and ( 10) ] was included in the computations. Since the atomic oxygen 
ion density-height profiles required for the calculations were not available, the 
0 + number density at each altitude, z, was computed from 
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(30) 


n[° + ] z = 


0. 95 x [n(E) ] 2 
z 


n(E) z + 


Yi n (N 2 ) 

Qfl 


Z ^ r 2 , . 

- + — X n(Oo) 

01 2 ^'Z 


in a manner similar to that of Wickwar [ 142] , where n(E) is the electron 

z 

number density at altitude z , n(N 2 ) and n(0 2 ) are the molecular nitrogen 

Z Z 

and oxygen number densities at altitude z , and the Greek symbols are the rate 
coefficients for the following reactions: 


0 + + N 2 -* NO + + N 

: Yl 

( 31 ) 

NO + + e — N + O0D) 

: cti 

(32) 

+ + 


(33) 

o + o 2 — o 2 + o 

' Y2 

0 2 + + 6 — 0 + 0* 

: oi 2 

(34) 


In this formulation O is assumed to be in a steady-state chemical 
equilibrium and it can never exceed 95 percent of the electron density, even at 
and above the F 2 peak. Neither of these assumptions is critical to the calculated 
production of 6300 A emission intensity. 

The number of *D oxygen atoms at altitude z from this production - 
mechanism is calculated by 

n[0( 1 D) ] = n[0 + ] [" n(0 2 ) x y 2 EFF1 + n(N 2 ) xy ( EFF2 1 

Z Z I Z Z 1 

(35) 

where EFF1 and EFF2 are the efficiencies of these reactions for producing *D 
oxygen atoms. 

The number density of 0(*D) s at altitude z is the total of these two 
production mechanisms. 
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C. Volume Emission Rate 


The volume emission rate at altitude z is computed from 
6. 9x 10 -3 x ntO^DjJ 

E 6300 - 9. 1 X 10“ 3 + x 'n(Nj) + q^x „(0 2 ) < ph0tOnS 0m ' S S ''> ' 

z z 

(36) 

where q t and q 2 are the rate coefficients for quenching by molecular nitrogen 
and oxygen, respectively. 


D. Integrated Emission Intensity 

Integrated emission intensity is calculated by using a trapezoidal integra- 
tion technique. 


VI. SOLAR FLUXES AND ABSORPTION CROSS SECTIONS 

Two separate sets of solar flux values in the Schumann-Runge continuum 
with associated O z absorption cross sections were used in this study. Table 13 
contains solar flux values from Friedman and Tousey of the Naval Research 
Laboratory as reported by Hinteregger, Hall, and Schmidtke [6]. The cross 
sections contained in this same table are from Metzger and Cook [7] . Table 
14 contains solar fluxes and absorption cross sections recommended by 
Ackerman [5] . 


VII. DISCUSSION 


On the night of January 6, 1972, conjugate point sundown occurred at 
approximately 2233 CST, assuming a screening height of 180 km; therefore, 
the airglow over the observation point (Huntsville, Ala. ) after this time could 
only be the result of dissociative recombination of the molecular ions, 0 2 + and 
NO + . On this night only the observations at 2303 and 2330 CST occurred after 
2233 CST. 
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TABLE 13. SCHUMANN-RUNGE CONTINUUM SOLAR FLUXES OF 
FRIEDMAN AND TOUSEY WITH ASSOCIATED O z ABSORPTION 
CROSS SECTIONS FROM METZGER AND COOK [6,7] 


Wave Length 

(A) 

Solar Fluxes 
(cm -2 s -1 ) 

Cross Sections 
(cm 2 ) 

1775 - 1725 

1.06 x 10 12 

1.70 x 10~ 19 

1725 - 1675 

7.00 x 10 11 

8.60 x io -19 

1675 - 1625 

4.20 x lo 11 

2.20 X 10' 18 

1625 - 1575 

2.60 x 10 11 

5.20 x 10" 18 

1575 - 1525 

1.33 x 10 11 

8.00 X 10 -18 

1525 - 1475 

7.30 x io 10 

1. 13 x 10 -17 

1475 - 1425 

3.70 x 10 10 

1.43 x 10 -17 

1425 - 1375 

1.83 x io 10 

1.43 x 10 -17 

1375 - 1325 

1 

1.77 X IO 10 

5.80 X 10 -18 


Airglow intensities for these two times were calculated using the set of 
equations discussed in detail in Section V with all three atmospheric models 
discussed in Section IV and the electron density-height profiles discussed in 
Section III. The neutral atmospheric constituent distributions are established 
by the derived exospheric temperatures. Once these components are established, 
the only inputs to the set of equations that affect the calculated intensities are 
the rate coefficients for equations (31) through (34) , the efficiencies of the 
0( J D) production mechanisms in equation (35) , and the quenching coefficients 
in equation (36) . 
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TABLE 14. SCHUMANN-RUNGE CONTINUUM SOLAR FLUXES AND 
ASSOCIATED O a ABSORPTION CROSS SECTIONS RECOMMENDED 

BY ACKERMAN [5] 


AA (A) 


Lya 1.215,67 
1.170 - 1.163 
1.176 - 1.170 
1.183 - 1.176 
1.190 - 1.183 
1.198 - 1.190 
1.205 - 1.198 
1.212 - 1.205 
1.220 - 1.212 
1.227 - 1.220 
1.235 - 1.227 
1.242 - 1.235 
1.250 - 1.242 
1.258 - 1.250 
1.266 - 1.258 
1.274 - 1.266 
1.282 - 1.274 
1.290 - 1.282 
1.299 - 1.290 
1.307 - 1.299 
1.316 - 1.307 
1.324 - 1.316 
1.333 - 1.324 
1.342 - 1.333 
1.351 - 1.342 
1.360 - 1.351 


q (cm -2 s -1 ) 


3.00 x 10 11 
1. 03 x 10? 


2.28 x io 9 


1.24 x IO 10 
1.99 x io 9 


a(0 2 ) (cm2) 


1.00 X IO -20 

2.00 x io -20 
1.25 x io -18 
2. 55 X 10’ 19 

3.00 x IO -20 

3.75 X 10’ 19 

4.45 x io' 18 

8.35 

6.00 x 10' 19 

2.35 
4.50 

3.35 

1.75 x 10“ 17 
8.95 x 10" 19 


4.20 


1.45 x io' 18 


2.30 x 10 18 


1.370 - 1.360 







TABLE 14. (Concluded) 


AA(A) 

q (cm -2 s -1 ) 

0 ( 02 ) (cm 2 ) 

1.379 - 1.370 

2.74 

1.23 X 10' 17 

1.389 - 1.397 

3.10 

1.32 

1.408 - 1.389 

7.60 

1.36 

1.428 - 1.408 

1. 01 x lo 10 

1.40 

1.449 - 1.428 

1.30 

1.48 

1.470 - 1.449 

1.82 

1.41 

1.492 - 1.470 

2.33 

1.29 

1.515 - 1.492 

2.66 

1.15 

1.538 - 1.515 

2.90 

9.91 x lo -18 

1.562 - 1.538 

3.60 

8.24 

1.587 - 1.562 

4.75 

6.58 

1.613 - 1.587 

6.40 

4.97 

1.639 - 1.613 

5.49 

3.45 

1.667 - 1.639 

1.19 x 10 11 

3.08 

1.695 - 1.667 

1.76 

1.23 

1.724 - 1.695 

2.32 

7.22 x 10 -19 

1.739 - 1.724 

1.44 

4.58 

1.754 - 1.739 

1.83 

2.74 


A. Reaction and Quenching Rate Coefficients 

It was assumed that the results of Dalgarno and Walker [11] were 
correct and that no O^D) atoms were formed from the dissociative recombina- 
tion of NO + . Since quenching of 0( 1 D) atoms by 0 2 is insignificant in compari- 
son with quenching by N 2 and its reaction rate coefficient appears to be reliably 
known (Table 7) , the value of 5 x 10” 11 cm 3 s -1 was used throughout the calcu- 
lations. With the exception of the N 2 quenching coefficient and the efficiency 
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of the 0 2 + e production mechanism, the remaining rate coefficients also 
appear to be reliably known; therefore, the following values were used throughout 
the calculations: 


+ . a 


« 1 

N + O 

«i = 

« 2 

O + O 

a 2 = 

72 

+ 


* 

o 2 +o 

72 = 

7l 

Hr 



NO + N 

7i = 


Q 2 ^ 


o 2 

- o 2 + o 

q 2 = 


OH = 4. 1 x 10 -7 (300/T) cm 3 s" 1 
a 2 = 1. 95X10" 7 (300/T) 0,7 cm 3 s -1 


«-i 


y t = 1 x 10 -12 cm 3 s _1 


i-ll ^w.3 „-l 


With these input parameters fixed, there are only two parameters that can be 
varied to obtain agreement between the calculated and observed intensities at 
2303 and 2330 CST. These are the N 2 quenching coefficient and the efficiency 
of the 0 2 + + e production mechanism. All possible combinations of these two 
rate coefficients were used with all three model atmospheres. In all of the 
succeeding figures the solid line connects the observed intensity values while 
the broken line connects the calculated intensity values. Figures 20 through 22 
show results using the Jacchia 1964 model atmosphere. 

In Figure 20 the molecular nitrogen quenching coefficient is 1 x 10 -1 ° 
cm 3 s _1 and the 0 2 + + e production mechanism efficiency is 1.00, which means 
that one 0( 1 D) atom is formed from every dissociative recombination. These 
are maximum values for both of these parameters, as reported in the literature. 
The calculated intensities at both 2303 and 2330 CST are approximately three 
times the observed intensities. Calculated intensities during the remainder of 
the period being studied are substantially larger than the observed intensities. 

In Figure 21 the quenching coefficient remains as lx 10 -1 ° cm 3 s' 1 while 
the dissociative recombination production efficiency is decreased to 0.50, which 
means that only one 0( 1 D) atom is formed from every two recombinations. The 
calculated intensities at 2303 and 2330 CST are still almost 50 percent larger 
than the observed intensities. A similar decrease is evident in the calculated 
intensities during the remainder of the period being analyzed, although they are 
still larger than the observed intensities. 
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SUNSET 


TIME (CST) 

Figure 20. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6r, 1972. 


In Figure 22 the quenching coefficient value remains at 1 x 10“ 10 cm 3 s -1 
while the efficiency of the dissociative recombination production mechanism is 
decreased to 0. 333, which means that only one O^D) atom is produced from 
every three recombinations. Agreement is excellent between calculated and 
observed intensities at both 2303 and 2330 CST. Once again, as expected, 
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Figure 21. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 


agreement has improved between calculated and observed intensities during the 
remainder of the period. For the Jacchia 1964 model atmosphere, then, the 
best combination of parameters is a quenching coefficient of 1 x 10 -1 ° cm 3 s -1 
and a dissociative recombination efficiency of 0. 333. This result tends to 
reinforce the assumption that the NO + + e mechanism does not produce any 
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Figure 22. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 


0( 1 D) atoms, since any additional O^D) atoms would require either a quenching 
coefficient larger than the value used, which is the maximum found in the litera- 
ture, or an 0 2 + + e production efficiency of less than 0. 333. This, however, 
contradicts the results of Wickwar [142] who found that a value of 0.50 provided 
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the best results. A larger quenching coefficient would not agree with the results 
obtained by Hernandez [113] from the Ionospheric Modification Experiment 


which yielded a value of 9 




lO -10 cm 3 s -1 based on a Jacchia 1964 model 


atmosphere formulation using an N 2 number density boundary value of 2. 0 x 10 11 
cm -3 instead of the 4. 0 x 10 11 cm r3 values used herein. 


Figures 23 and 24 show results using the Jacchia 1971 model atmosphere. 
In Figure 23 the quenching coefficient is 5 x 10 -11 cm 3 s - ' 1 , and the dissociative 
recombination production efficiency is 0. 50. The calculated intensities at both 
2303 and 2330 CST are approximately 60 percent larger than the observed 
emission intensities. Calculated intensities are greater than the observed 
intensities throughout the observation period. In Figure 24 the quenching 
coefficient is increased to 1 x 10 _1 ° cm 3 s" 1 . Agreement between calculated 
and observed intensities at both 2303 and 2330 CST is perfect. Calculated 
intensities are greater than observed intensities during that portion of the period 
when sunlight contributes to the emission intensity, prior to approximately 
1805 CST, and less than observed intensities from approximately 1805 to 2303 
CST. Any other combination of values of these parameters would result in the 
use of either a value of the quenching coefficient larger than the maximum value 
quoted in the literature or a value of the dissociative recombination production 
mechanism efficiency of less than 0. 50, the value currently quoted in the 
literature as being most acceptable. 

Figures 25 through 27 show the results of these combinations using the 
120 km boundary conditions recommended by Reber 3 with the Jacchia 1964 model 
formulation and a constant temperature shape parameter of 0.01875, also 
recommended by Reber. These values were established by Reber in conjunction 
with his efforts to establish a new U.S. Standard Atmosphere; however, the 
resulting number density-height profiles he established are significantly 
different from those calculated herein because he used vertical flux and molecu- 
lar diffusion terms that effectively changed the models from static to pseudo- 
dynamic. These fictitious flux terms, which were used to ensure compatibility 
between the model and observational data, increased the molecular nitrogen 
number densities while decreasing the molecular oxygen number densities at 
all altitudes. Both of these modifications would result in a decrease in the 
calculated intensity of the emission line. 

In Figure 25 the quenching coefficient is 5 x 10 -11 cm 3 s -1 , while the 
dissociative recombination production mechanism efficiency is 0.50. Calcu- 
lated intensities are almost 80 percent larger than the intensities observed at 
both 2303 and 2330 CST. Calculated intensities are larger than the observed 
intensity value throughout the period. 


3. Ibid. 
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Figure 23. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 


In Figure 26 the quenching coefficient has been increased to 1 x 10 10 
cm 3 s -1 without changing the efficiency of the dissociative recombination pro- 
duction mechanism. Agreement between calculated and observed emission 
intensities has been improved throughout the period; however, the calculated 
int ens ities are still larger than those observed at both 2303 and 2330 CST. 
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SUNSET TIME (CST) 

Figure 24. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 


In Figure 27 the quenching coefficient remains at 1 x 10 10 cm 3 s -1 , while 
the production mechanism efficiency has been reduced to 0. 434. Agreement 
between calculated and observed intensities at both 2303 and 2330 CST is 
excellent. That is the best combination of values for this model atmosphere; 
however, if the fictitious flux and diffusion terms recommended by Reber had 
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SUNSET TIME ( CST ) 

Figure 25. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 


been included so that his number density-height profiles would have been pro- 
duced by the model, a quenching coefficient of approximately 9 x 10 -11 cm 3 s -1 
and a production mechanism efficiency of approximately 0.60 would probably 
have been the best combination of values. 
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Figure 26. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 


B. Solar Radiation Fluxes and Associated Absorption 
Cross Section Effects 

Once agreement between calculated and observed intensities resulting 
only from the dissociative recombination process has been attained, it seems 
reasonable to assume that the preferred values of the quenching rate coefficient 
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Figure 27. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 


and the 0 2 + + e production mechanism efficiency so established for each atmos- 
pheric model should remain constant during the remainder of the comparisons 
in which other parameters make an important contribution. In this section the 
Schumann-Runge continuum fluxes and their associated absorption cross sections 
are varied to study their effects on the calculated intensities. 
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Figures 28 and 29 are comparisons of observed intensities with intensities 
calculated using the Jacchia 1964 model atmosphere, a quenching rate coefficient 
of 1 x 10 -I ° cm 3 s" 1 , an O z + + e production mechanism efficiency of 0. 333, and 
Hinteregger et al. [6] reported fluxes with the Metzger and Cook [7] cross 
sections (Table 13) and the Ackerman [5] recommended fluxes and cross sections 
(Table 14), respectively. 



Figure 28. 6300 A emission line. intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 
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Figure 29. 6300 A emission line intensities derived from Fabry- Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on J anuary 6, 1972. 


In Figure 28 the agreement of the calculated and observed intensities at 
2303 and 2330 CST deteriorates immediately and becomes progressively worse 
as one gets nearer the time of sunset. Between sunset and approximately 1811 
CST, the time period when the dissociative photoexcitation process contributes 
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significantly to the overall intensity of the 6300 A emission, the calculated 
intensities are significantly larger than the observed intensities. From 1811 
CST to approximately 2300 CST, the calculated intensities are smaller than the 
observations with the difference between the two decreasing with time after 
sunset. 


In Figure 29 the agreement of the calculated and observed intensities at 
2303 and 2330 CST once again becomes progressively worse the nearer one gets 
to the time of sunset; however, the major difference between these results and 
those shown in Figure 28 is that the calculated intensities are always smaller 
than the observed intensities. 

From a comparison of these two analyses, it appears that either the 
Hinteregger et al. reported fluxes are not representative of the incoming solar 
flux on the night of January 6, 1972, or the associated Metzger and Cook cross 
sections are incorrect. It also appears that if the Ackerman recommended 
fluxes and absorption cross sections are indeed representative of the solar con- 
ditions, there is another source of 0(*D) atoms on this night. The magnitude 
of this contribution is shown as a function of time in Figure 30. This is com- 
parable to the results reported by Wickwar [36] and attributed to direct impact 
excitation of atomic oxygen by conjugate point photoelectrons. His results and 
those here show that the enhancement continues until at least conjugate point 
sundown which occurred on January 6, 1972, at approximately 2233 CST assuming 
a screening height of 180 km. At this time the conjugate point solar zenith angle 
was approximately 103°, which is less than the approximate 105° conjugate point 
solar zenith angle at which Wickwar [36] reported a cessation of the enhancement 
over Arecibo, P.R. The time at which the conjugate point (CP) photoelectron 
contribution ceased is in excellent agreement with Figure 5 of Reference 143 
which shows that during a period of low solar activity, CP photoelectrons should 
not be present over the Huntsville area after 2300 CST in January. The magni- 
tude of this conjugate point photoelectron impact excitation contribution is much 
larger than the approximate 20 R contribution over Blue Hill Observatory, Mass. , 
reported by Noxon and Johanson [ 133] and somewhat larger than that reported 
by Wickwar [36]. 

Figures 31 and 32 are comparisons of observed intensities with calculated 
intensities using the Jacchia 1971 model atmosphere, a quenching rate coefficient 
of 1 x 10 -10 cm 3 s -1 , and a dissociative recombination production mechanism 
efficiency of 0. 50 with the Hinteregger et al. [6] and Metzger and Cook [7] 
combination of fluxes and cross sections and the Ackerman [5] fluxes and cross 
sections, respectively. The only differences in these two calculations occur 
prior to 1805 CST where the intensities calculated with the Hinteregger et al. 
and Metzger and Cook data combination (Fig. 31) are approximately 80 percent 
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Figure 30. Difference between observed and calculated 6300 A emission line 
intensities as a function of time after sunset, attributed to direct impact 


excitation of atomic oxygen by conjugate point photoelectrons. 


larger than the observed intensities, while those calculated with the Ackerman 
data are about 20 percent less than the observed. If the Jacchia 1971 model is 
truly representative of the atmosphere on this night, the Hinteregger et al. and 
Metzger and Cook combination is not a representative set of solar fluxes and 
absorption cross sections. On the other hand, if the Jacchia 1971 model and 
the Ackerman data are both representative, there must be a third source of 
0( 1 D) atoms. Again, this could be the direct impact excitation of oxygen atoms 
by conjugate point photoelectrons with the contribution at Huntsville continuing 
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Figure 31. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 


until conjugate point sundown at approximately 2233 CST. Figure 33 shows the 
magnitude of this third contribution as a function of time. As with the results 
obtained using the Jacchia 1964 model, the contribution ceases approximately 
at the same cutoff time observed over Arecibo, P. R. , as reported by Wickwar 
[36]. 
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Figure 32. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 


Figures 34 and 35 show the same comparison using the Jacchia 1964 
model formulation with the boundaiy conditions and temperature shape parameter 
proposed by Keber without the vertical diffusion and flux corrections. Although 
these results can be discussed only qualitatively, they are identical to those 
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Figure 33. Difference between observed and calculated 6300 A emission 
line intensities as a function of time after sunset, attributed to direct 
impact excitation of atomic oxygen by conjugate point photoelectrons. 


using both the Jacchia 1964 and 1971 model atmospheres. If both this model and 
the Ackerman data are representative, then a third source of 0( J D) atoms of 
approximately the same magnitude is present with the variation with time as 
shown in Figure 36. 
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Figure 34. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 


C. Neutral Gas Temperature Effects 

Estimated errors in temperatures derived from the F-P observations 
range from approximately ±85 K near sunset to approximately ±195 K at 
midnight, as shown in Figure 18. If the actual temperatures are higher than 
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Figure 35. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 


the ones derived from the Fabry-Perot interferometer measurements that were 
used in theoretically calculating the intensities used in this analysis, then: 

1. The N 2 quenching rate coefficient of 1 x 10 _1 ° cm 3 s -1 , determined 
as a result of this analysis, is too small by approximately 9 to 10 percent, or 
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Figure 36. Difference between observed and calculated 6300 A emission line 
intensities as a function of time after sunset, attributed to direct impact 
excitation of atomic oxygen by conjugate point photoelectrons. 

2. The dissociative recombination production mechanism efficiency of 
0. 50 is too large by approximately 9 to 10 percent, or 

3. Both of these results are in error by a lesser amount such that the 
combination is in error by approximately 9 to 10 percent. 

Figures 37 and 38 are plots of constant N 2 number density and volume 
emission rate, respectively, as a function of solar zenith angle for January 
6-7, 1972, for each time during the observational period when simultaneous 
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observations of the electron density-height profile, the integrated intensity and 
an exospheric temperature from the Fabry-Perot interferometer, and integrated 
intensity from the turret photometer were available. The solid line in each of 
the figures shows the computational results using the actual temperatures derived 
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Figure 38. Contours of constant volume emission rate. 


from the F-P measurements, while the dotted lines show the results using 
temperatures increased by 5Q K. Figure 37 shows that the model of the atmos- 
phere expands upward with increased exospheric temperature, while Figure 38 
shows that the contours of constant volume emission rate are also displaced 
upward while the peak volume emission rate values are increased approximately 
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2 to 3 percent. The calculated integrated intensities are increased 9 to 10 
percent. Figure 39 shows calculated and observed intensities as a function of 
A im? with the same inputs to the calculations, while Figure 40 shows the third 
source contribution under these conditions. Exospheric temperatures derived 
from Fabry-Perot interferometer measurements are used in these calculations 



Figure 39. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 
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Figure 40. Difference between observed and calculated 6300 A emission line 
intensities as a function of time after sunset, attributed to direct impact 
excitation of atomic oxygen by conjugate point photoelectrons. 


with the Jacchia 1971 model atmosphere, a quenching rate coefficient of 1 x 10~ 10 
cm 3 s -1 , a dissociative recombination production mechanism efficiency of 0. 50, 
and the Ackerman [ 5] recommended values for the Schumann-Runge continuum 
with associated absorption cross sections. In these figures the altitude of the 
peak of the volume emission rate, I, can be considered to be the same as the 
altitude of the derived temperature. The height of this peak is variable; how- 
ever, in all instances it is between 220 and 255 km. In the Jacchia 1971 model 
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atmosphere this is well below the altitude where the atmosphere becomes iso- 
thermal, for all practical purposes; therefore, the derived temperatures should 
be lower by 15 to 25 K than the true exospheric temperatures. 

Figure 41 shows the observed and calculated intensities on January 6, 
1972, where the calculations are based on the same inputs as those used for the 
results shown in Figure 39 with the exception that the exospheric temperature 
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Figure 41. 6300 A emission line intensities derived from Fabry-Perot 
interferometer observations and theoretically calculated intensities 
as a function of time after sunset on January 6, 1972. 
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values are all increased 50 K. The third source contributions under these con- 
ditions are shown in Figure 42. Since the contributions at 2303 and 2330 CST 
are both negative, the dissociative recombination process efficiency would have 
to be reduced and, once this is done, the results would be identical to those 
shown on Figure 40. Similar calculations have been made for the other models; 
however, since the results are almost identical, they are not shown here. 
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Figure 42. Difference between observed and calculated 6300 A emission line 
intensities as a function of time after sunset, attributed to direct impact 
excitation of atomic oxygen by conjugate point photoelectrons. 


81 



D. Electron Density-True Height Profile Effects 


There are several problems associated with the true height analysis of 
ionograms which, although they are outside the scope of this analysis, affect the 
results. Problems associated with the reduction of daytime ionograms obviously 
do not affect the results of this analysis; however, the start problem discussed 
in detail in ionospheric literature, as well as all the scale problems, can sig- 
nificantly affect the results. 

Scaling problems arise from multiple traces resulting from multiple 
reflections, some of which, or perhaps all, are from off- vertical reflections, 
proper selection of points to be used in the program, correct tracing of the 
negative on the scaling form, and correct identification of the ordinary and 
extraordinary ray signals. From the ionograms created by the standard C-4 
ionosonde used in this study, it is impossible to discriminate between vertical 
and off-vertical returns if there are multiple traces on the film. Fortunately, 
for this analysis, enough ionograms were available to preclude the use of any 
with multiple traces; however, there is no way of knowing whether the single 
traces on the ionograms are from vertical incidence reflections or not. They 
have all been assumed to be from vertical incidence reflections. 

As previously stated in Section in, the start problem was investigated 
by using several different programs in conjunction with several known electron 
density-true height profiles until the data reduction program output the known 
electron density-true height profiles, to the accuracy required, from which the 
input data were originally selected. This technique should minimize the 
possibility of large errors in the electron density-true height profiles; however, 
if errors do occur, the analysis shows that the reduction technique tends to 
underestimate the height of the F 2 maximum. These errors in height can be on 
the order of 10 to 20 km with the mean errors expected to be closer to 10 km. 

If the actual height of the F 2 maximum is approximately 10 km higher than the 
height calculated by the true height reduction technique, that portion of the 
intensity due to dissociative recombination computed using the calculated height 
will be approximately 20 percent greater than it actually is. This will lead to 
overestimating the value of the quenching coefficient and/or the molecular 
nitrogen number density and/or underestimating the dissociative recombination 
production mechanism efficiency and/or the molecular oxygen number density 
or various combinations of these. 

Accurate electron density-true height profiles, especially near the F 2 
maximum, are a necessity in 6300 A emission line intensity analyses. 
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E. General 


Figure 43 is a composite of the observational and theoretical data for 
the night of January 6, 1972. From Figures 43(a) and (b) it is readily apparent 
that each time the height of the F 2 maximum changed, there was a corresponding 
change in the temperature derived from the Fabry-Perot interferometer. The 
second line on Figure 43(a) shows that the altitude of the calculated maximum 
volume emission rate also varied in agreement with the two measured quantities; 
therefore, the temperature of a higher altitude layer of the atmosphere is being 
measured. Correspondence between changes in the altitude of the F 2 maximum 
and the height of the peak value of the volume emission rate is not as good 
during the early portion as it is during the latter portion of the twilight period. 
This is because there is a substantial contribution from the dissociative photo- 
excitation process in addition to the dissociative recombination process. This 
supports the conditions of Roble, Hays, Nagy [124] that the 0( 1 D) atoms are 
in thermal equilibrium with the neutral atmosphere and 6300 A emission can be 
used to monitor the exospheric temperature. 

Figure 43(c) shows the electron density at the F 2 maximum, while Figure 
43(d) shows a combination of observed and calculated intensities as a function 
of solar zenith angle. The top curve in Figure 43(d) represents the total 
observed intensity on the night of January 6, 1972. The bottom curve in this 
figure connects the intensity values attributed from the theoretical calculations 
at each observation time to the dissociative recombination of molecular oxygen 
ions. 


On the left side of Figure 43 the difference between the short horizontal 
line and the observed intensity curve is attributed, from the calculations, to the 
dissociative photoexcitation of molecular oxygen by solar fluxes in the 
Schumann-Runge continuum. The remainder of the observed intensity, the 
difference between the bottom curve and the curve composed of the short hori- 
zontal lines at a value of 183 R and the total intensity curve between solar 
zenith angles of 104° and 162° , is believed to be due to direct impact excitation 
of atomic oxygen by conjugate point photoelectrons. As expected, there appears 
to be an excellent correspondence between the electron density curve and the 
curve that represents the intensity attributable, from the calculations, to the 
dissociative recombination process. This is shown in more detail in Figure 44 
where the O *s show the observed intensity as a function of the electron density 
at the F 2 peak. There appears to be an excellent correspondence between the 
two parameters with the exception of those three points relating the higher 
intensities to the larger electron densities, and the deviation here is due to the 
contribution to the intensity from the dissociative photoexcitation process. If 


83 


ALT 260- 
(km) 



• F2 MAXIMUM 


MAXIMUM 

VOLUME EMISSION 
RATE 



EXOSPHERIC 
TEMPERATURES 
DERIVEO FftOM 
OBSERVATION 


F 2 MAXIMUM 
ELECTRON DENSITY 


CONJUGATE POINT SUNSET 
ASSUMING 180 km 
SCREENING ALTITUDE 


PORTION OF OBSERVED INTENSITY 
CALCULATED TO BE DUE TO: 

V//77//2 PHOTO DISSOCIATION 
miimill IMPACT EXCITATION 
Enzzza DISSOCIATIVE RECOMBINATION 


OBSERVED INTENSITY 


104 112 


120 128 


SOLAR 2ENITH ANGLE 


Figure 43. Composite of observational and theoretical data as a function 
of local solar zenith angle on January 6, 1972. 



INTENSITY (R) 

Figure 44. Observed and calculated 6300 A intensities versus the electron 

number density at the F 2 peak. 

this intensity was adequately calculated and its contribution removed, it would 
seem as if the resulting intensities should fall on the straight line defined by 
the other observed intensities and electron densities. The squares, which 
denote the results of subtracting the intensity calculated to be due to dissociative 
photoexcitation from the observed intensity, do not indicate the same correspond- 
ence between the electron density and 6300 A intensity as the other observations; 
therefore, it appears as if there is a possibility that either the dissociative 
photoexcitation process is improperly modeled or there is a change in the elec- 
tron density-dissociative recombination intensity relationship when the sun is 
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contributing to the emission line intensity. These results are identical to those 
shown in Figure 43 (d) , and both seem to indicate that the dissociative photo- 
excitation and dissociative recombination processes combine in some manner 
to effectively control the effect of conjugate point photoelectrons on the local 
6300 A OI emission. On Figure 44 the circles show the intensities, calculated 
on the basis of the dissociative recombination process only, as a function of the 
electron density at the F 2 peak. Correspondence between these two parameters 
is excellent, as ejected. The equation of a straight line fit to these data is 

1(B) = (y) (N e -2.50) x 10-4 . 

When the calculated dissociative photoexcitation process contribution was sub- 
tracted from each observed intensity when the sunlight was making a contribu- 
tion, a value of 183 R remained in each instance. From the previous equation 
the intensity should be 183 R when the electron density is 6.672 x 10 5 cm -3 . 
Perhaps this is an indication of how the conjugate point effects are controlled 
by the local ionosphere, but verification of that is beyond the scope of this 
analysis and, undoubtedly, would require much additional observational data. 

Figure 45 compares the 6300 A total emission intensity over Huntsville, 
Ala. , on the night of January 6, 1972, with the mean winter twilight 6300 A 
intensities over the Blue Hill Observatory, Mass. , during the 1968-1970 time 
period. During the early portion of the twilight period, the Blue Hill Observatory 
intensities are about 50 percent greater than those over Huntsville. During the 
latter part of the twilight period, the difference between the two sets of intensity 
observations increases to a factor of 3. 5 to 4. Figure 46 provides a partial 
explanation for this difference in observed intensities in that the mean winter 
time densities at the F 2 peak over Blue Hill for the 1968-1970 time period are 
between 1. 5 and 2. 5 times as large as those over Huntsville on January 6, 1972. 


VIII. CONCLUSIONS 


A. Atmospheric Model 

The results of this analysis indicate that the Jacchia 1971 model atmos- 
phere is more representative of the molecular oxygen and nitrogen concentra- 
tions over Huntsville, Ala. , on the night of January 6, 1972, than other existing 
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Figure 45. Comparison of 6300 A emission intensities over Huntsville, Ala, 
and Blue Hill Observatory, Mass. , as a function of solar zenith angle. 
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Figure 46. Comparison of F 2 peak electron number densities over 
Huntsville, Ala. , and Blue Hill Observatory, Mass. 



models. All other models require either dissociative recombination production 
mechanism efficiencies lower than the currently preferred value of 0. 50 or rate 
coefficients for the quenching of 0( 1 D) by molecular nitrogen higher than the 
highest value currently accepted by the scientific community, 1 x 10" 10 cm 3 s -1 . 
The imposition of these limits on theoretical calculations of the 6300 A emission 
shows that the Jacchia 1971 model consists of a lower limit on the molecular 
nitrogen number density and an upper limit on the molecular oxygen number 
density. Current planned revisions to the Jacchia 1971 model atmosphere 
include an increase in the molecular nitrogen number density and a decrease in 
the molecular oxygen density. These changes could readily be accommodated 
within the rate coefficient limits discussed above. 

The 6300 A emission can provide data concerning the N 2 and 0 2 densities 
for use in the refinement and development of atmospheric models; however, 
additional emission lines must be monitored simultaneously to obtain the data 
required for the future refinement of the models at orbital altitudes where the 
principal constituent is atomic oxygen. 


B. Photochemical Rate Coefficients 

Results show that the following values constitute an internally consistent 
set that is capable of describing the post-sunset twilight behavior of the 6300 A 
emission intensity: 


+ + 

o + o 2 — *■ o 2 + o 

0 2 + + e — O* + O 


0 + + N 2 — NO + + N 
NO + + e — N + O* 

O^D) + O z — O + O z 
0( 1 D) +N 2 ^0 + N 2 


2 x 10 -11 cm 3 s -1 


1. 95 x 10" 7 (300/T) 0,7 cm 3 s _i with an 
efficiency of producing 0( I D) atoms of 


0.50 


+0.32 

- 0.12 


1 x 10 -12 cm 3 s -1 

4.5 x 10“ 7 (300/T) cm 3 s” 1 with an 
efficiency of producing O^D) atoms 
of 0.00 

5 x 10 -11 cm 3 s -1 



cm 3 s -1 
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The first three of these are upper limits since use of larger values would require 
that the upper limit on the N 2 quenching rate coefficient be exceeded or the 
dissociative recombination mechanism efficiency be less than 0. 50, while the 
last two values are those currently most acceptable. 


C. Solar Flux and Absorption Cross Sections 

The conclusions reached above rule out the possibility of the Friedman 
and Tousey solar flux data as reported by Hinteregger et al. [6] as being 
representative of the solar conditions over Huntsville, Ala., on January 6, 1972, 
if the data are to be used in combination with the Metzger and Cook [7] absorp- 
tion cross sections. The Ackerman [5] recommended solar fluxes and absorp- 
tion cross sections are acceptable; however, additional observational data, such 
as fluxes and energy spectra of conjugate point photoelectrons, are required to 
positively establish the exact combination of solar flux and absorption cross 
section at any single point in time. For example, the Schumann-Runge contribu- 
tion at 1726 CST could be increased from 67 R to 109 R and the agreement 
between the observed and calculated intensities would be perfect without the 
requirement for a large conjugate point photoelectron contribution. Therefore, 
although the Schumann-Runge fluxes and associated absorption cross sections 
listed in Table 14 are acceptable, an upper limit could be established as that 
combination of fluxes and absorption cross sections which would produce 109 R 
of 6300 A OI intensity at 1726 CST. However, at 1726 CST the F 2 maximum 
electron density is decreasing rapidly with time and it is reasonable to assume 
that at least a comparable rate of decrease existed between sunset at 1703 CST 
and 1726 CST. Therefore, it is believed that this upper limit would not be 
realistic, since it appears that its use at smaller solar zenith angles would 
result in the calculated intensities becoming larger than the observed intensities. 

Since no observational data are available prior to a local solar zenith 
angle of 97.5°, this analysis cannot provide any additional information about the 
conclusion of Noxon and Johanson [133] that the local photoelectron direct 
impact excitation process must be included for solar zenith angles between 90° 
and 95°, since the dissociative photo excitation process alone cannot account for 
the difference between observed intensities and those calculated from the 
dissociative recombination process alone. In the winter the effects of CP 
photoelectrons further complicate the analysis. 
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D. Conjugate Point Photoelectrons 

There is a sizeable contribution to the 6300 A emission line intensity 
from conjugate point photoelectrons, at least during the period from a local 
solar zenith angle of approximately 97. 5° to a conjugate point solar zenith angle 
of approximately 103°. These limits are set in this analysis by the earliest 
observation time and conjugate point sunset, assuming a screening height of 
180 km. Based on the calculations in this analysis, the contribution of the 
combined dissociative recombination and conjugate point photoelectron direct 
impact excitation processes to the 6300 A emission line intensity was approxi- 
mately 183 R during the period from a local solar zenith angle of approximately 
97.5° to a local solar zenith angle of approximately 103°. If these observations 
and calculations can be extrapolated to earlier times in the twilight period, 
perhaps the results imply that conjugate point photoelectrons do not contribute 
to the 6300 A intensity as long as the local thermal electron densities are capable 
of producing 183 R through the dissociative recombination process. 


E. Atmospheric Motions 

Variations in the heights of the F 2 peak in the ionosphere were not 
associated with variations in the observed intensity; however, they were highly 
correlated with variations in the temperature derived from the Doppler broaden- 
ing of the 6300 A emission line. This implies that both the neutral and charged 
particle atmospheres were moving together and that observations of intensity 
enhancement by other investigators were more probably associated either with 
movements of the ionosphere that were not in phase with the neutral atmospheric 
movement or with travelling patches with increased electron densities. 

Even though the variations in the temperature observations were highly 
correlated with the height of the F 2 peak, the magnitudes of the variations are 
such that the changes cannot be accounted for solely on the basis of observing 
another altitude region in a quiescent atmosphere. Changes of 20 to 40 km in 
the altitude of the F 2 peak are associated with temperature changes on the order 
of 100°K. At the heights at which these observations apply, temperature-height 
profiles in current models are nearly isothermal and changes of 10°to 5Cf K 
would be expected over these height intervals. Therefore, either the 
temperature-height profiles in the models are in error or the observations are 
being made at different altitudes in air masses with different exospheric temp- 
eratures. Wind measurements from the Doppler shift of the 6300 A emission 
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line could be used to determine whether these temperature variations should be 
attributed to vertical and/or horizontal motions or to differences between the 
true and the model temperature-height profiles. No detailed measurements of 
the variations in the outside pressure, the reference pressure for the etalon 
chamber, are available for the time period being analyzed; therefore, no data 
are available that can be used to derive vertical wind speeds which might provide 
some insight into this problem. 


F. General 


The total 6300 A emission intensity over Huntsville, Ala. , appears to 
be much lower than that reported by other investigators [133, 142, 17] ; however, 
the conjugate point photoelectron contribution appears to be much larger during 
this declining portion of the 20th solar cycle. The lower integrated intensity 
over Huntsville can be attributed to the lower overall electron densities over 
Huntsville due to (1) its geographic location with respect to the geomagnetic 
equator which places it in an area of relatively low electron densities, (2) the 
fact that this analysis covered a twilight period in January when the annual 
variation would predict a minimum in the electron density, and (3) a time with 
respect to the solar cycle when the electron densities would be expected to be 
generally lower than those present for the time periods covered in the analyses 
referenced above. All of these factors would contribute to the increased CP 
photoelectron contribution observed over Huntsville if the local ionosphere 
controls the CP photoelectron input. 


IX. RECOMMENDATIONS FOR FUTURE RESEARCH 


The 6300 A emission line intensity and line shape are excellent tools for 
studying the Earth’s upper atmosphere, its composition, temperature, and 
dynamics; however, results are extremely sensitive to the electron density- 
height profiles used in the theoretical calculations required to substantiate and 
interpret the observational data. The phase path sounder at Huntsville, Ala. , 
which is currently being modified, will soon be operational and capable of 
producing electron density -true height profiles with an accuracy of ±0. 1 km. 

The increased accuracy of these data will allow a better determination 
of molecular oxygen and nitrogen density-height profiles, if independent experi- 
ments such as the Ionospheric Modification Experiment can provide reliable 
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results for the molecular nitrogen and/or oxygen quenching rate coefficients. 
Therefore, as soon as this sounder is available, the study should be reaccomp- 
lished to confirm these results. 

The phase path sounder and a CW Doppler array currently being installed 
will be capable of providing additional data concerning the horizontal and vertical 
motions of the ionosphere. The etalon chamber pressure should be referenced 
to a baratron so that line of sight motions of the 0( 1 D) atoms could be measured. 
The motions of both the charged particles and the 0( 1 D) atoms could then be 
used in conjunction with changes in the 6300 A emission intensity and the tem- 
perature of the neutral gas to study the effects of these motions on the emission 
intensity and the neutral gas temperature. 

Observational data for other periods should be subjected to this same 
analysis to ascertain the magnitude of the CP contribution as well as any possible 
seasonal variations in onset time and/or magnitude. Seasonal variations due to 
changes in electron densities should also be studied in conjunction with the 
observed variations in the emission line intensity. 

In view of the recent results of Hernandez 4 concerning the possible con- 
tamination of 6300 A OI emission observations by the ( 9-3) hydroxyl band, 
extreme care should be taken to ensure that the data used in any analysis 
consist solely of 6300 A photons from the atomic oxygen transition. Depending 
on its position with respect to the 6300 A OI line shape profile, the (9-3) OH 
band emission can produce changes that imply: 

1. Either lower or higher than actual temperatures for the 1 D oxygen 

atoms. 

2. More or less than actual atomic oxygen red line emission intensity. 

3. Mass motions of the *D oxygen atoms either toward or away from the 
observer. 

The results of Sivjee [144] are quite significant with respect to the 
6300 A OI emission intensities and derived exospheric temperatures obtained 
with the Fabry-Perot interferometer in its present configuration. The low and 
rapidly fluctuating temperatures derived from the emission line profiles on 
other nights might be due to changes in the intensity and temperature of the 
9-3 band OH emission (see Appendix D) . The F-P configuration should be 


4. Hernandez, op. cit. 
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changed so that emission lines of this band system do not fall at the peak of the 
6300 A OI emission line profile. In addition, other measurements, perhaps with 
photometers containing filters that would provide data on the P t ( 3) line of the 
9-3 OH band, should be made simultaneously so that the OH emission could be 
separated from the 6300 A OI emission. This combination of observational data 
could then be used to study the dynamics of the D -region of the atmosphere. 

This might resolve questions concerning the origin of the low temperatures 
occasionally derived from 6300 A OI emission line profiles. 
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APPENDIX A 


HISTORICAL BACKGROUND 


Newcomb [145] attributed the light of the night sky to the flux of radiation 
from stars that were too faint to be seen individually. Atmospheric light was 
studied intensively at the Lowell Observatory, Flagstaff, Arizona, from 1915 
onward, and the two atomic oxygen redlines, 6300 and 6364 A , were discovered 
there [146, 147], Seasonal, day-to-day, and daily variations of the airglow 
were gradually discovered [148-152] . Chapman [139] first proposed the theory 
that the airglow was caused by emission of photons associated with the recom- 
bination of ionized and, more particularly, dissociated particles produced by 
the absorption of sunlight, in the ultraviolet portion of the spectrum, during the 
day. Dufay [ 153] showed that the luminance of the night sky must come from a 
combination of the zodiacal light and the self-luminescence of the atmosphere, 
not from faint stars. 

Excitement of the 6300 and 6364 A atomic oxygen lines during twilight 
was first discovered by Garrigue [ 154, 155] . His investigations showed that 
after twilight the 6300 A emission diminished in intensity from 70 to 5 ( arbitrary 
units) at sunset at approximately 10° above the horizon and from 20 to 5 in the 
north at the same elevation angle. He concluded that the radiation was emitted 
by excitation resulting from solar radiation because of the progressive decrease 
after twilight. 

Cabannes and Garrigue [156] and Currie and Edwards [157] were the 
first to document and discuss the post-twilight decrease in the intensity of the 
6300 A emission. Their analysis of observations in December 1935 and January 
1936 at Pic du Midi and the spring of 1936 at Montpelier showed that the 6300 A 
emission intensity was affected by rapid perturbations, that it came from an 
altitude between 105 to 110 km, and that approximately 3 h after sunset it 
reached and thereafter remained fixed at a value approximately one-tenth the 
value at twilight with the intensity at the zenith two to three times weaker than 
at the horizon. These phenomena could not be explained by the action of sunlight 
even though they originated in the same atmospheric layers. The first explana- 
tions of the variation in intensity with the time based on the resonance scattering 
of sunlight were untenable; however, the photochemical dissociation of molecular 
oxygen by ultraviolet (X < 1750 A ) could not be truly effective at very high 
altitudes in view of the rarity of molecular oxygen [ 158] and the slow decrease 
after sunset required the presence of a very large number of oxygen atoms at 
very great heights, a very large scale height. 
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The spectrographic studies of Elvey and Farnsworth [159] and Elvey 
[160] showed that the slow decrease of the intensity of the 6300 A emission 
extended well into the night. The studies of Elvey and Farnsworth [159] also 
placed the height of the emitting layer at approximately 500 km with the twilight 
enhancement due to illumination of the atmosphere by the Sim, as postulated by 
Cabannes and Garrigue [156] . However, they also concluded that this effect 
could not account for all of the enhancement since it lasted for approximately 4 h 
after sundown. Resonance seemed unlikely because the emission was from a 
forbidden transition. Their study showed a difference between the morning and 
evening twilights; therefore, they concluded that the enhancement was due to 
some phenomenon of the upper atmosphere independent of the direct illumination 
of the Sun. From their pre-dawn enhancement analysis they concluded that 
there was a much slower decrease in the number of oxygen atoms with height 
than previous theoretical calculations had indicated. A gradual increase in 
intensity long before sunrise was also discovered at this time. 

Bates [161] postulated that the 6300 A emission originated in the F- 
region of the ionosphere and suggested that dissociative recombination of ions 
and electrons could be the source. Dufay and Tcheng Mao-Lin [162] confirmed 
both of these post- and pre-twilight enhancements that were investigated in more 
detail photo electric ally by Barbier [163] and photographically by Berthier [164] 
and Robley [165], Berthier’s [164] observations showed a very marked 
assymmetry between the morning and evening twilights. In the evening there 
was a very rapid decrease until the Sun was approximately 14° below the horizon 
and then a much slower decrease. According to his analysis the peak of the 
intensity occurred at an altitude between 90 and 100 km. 

It had seemed quite likely [ 166] that a large part of the twilight 6300 A 
emission arose from resonance scattering of sunlight; however, Chamberlain 
[ 15] showed that a comparison of observations with theory suggested that 
virtually none of the emission came from resonance scattering, a conclusion 
reached by Elvey and Farnsworth almost 15 years earlier. Chamberlain also 
showed that dissociation of molecular oxygen in the Schumann-Runge continuum 
contributed, but did not give by itself, either the observed intensity or its 
variations with time. He then advanced the hypothesis that the main excitation 
of the 6300 A twilight emission arose from dissociative recombination of 0 2 + 
following suggestions of Bates and Massey [61] , Bates [166] , and Barbier [167] ; 
however, he assumed that molecular oxygen was the principal deactivating 
particle. His study supported the results of Dufay and Tcheng Mao-Lin [168] 
and Barbier [ 14] that the emission height was approximately 250 to 300 km. 
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The rocket results of Heppner, Stolarik, and Meredith [ 169] and 
Heppner and Meredith [170] placed the emission height above 163 km, while the 
results of Zipf and Heath [171] raised the emission height to above 208 km. 

The rocket-borne measurements of Gulledge and Packer [ 172] and the satellite 
measurements of Reed and Blamont [173] gave further evidence that the maxi- 
mum of the 6300 A emission is at an altitude of approximately 250 to 300 km with 
the emitting layer on the order of 80 to 100 km thick. 

The ionospheric analyses of Ratcliffe, Schmerling, Setty, and Thomas 
[174] and Bates [175, 176] presented further evidence of the importance of 
dissociative recombination. St. Amand [177] and Barbier [167] showed a 
connection between the F-layer and the redlines. Huruhata, Nakamura, Tanabe, 
and Tohmatsu [178] concluded from their analysis of rocket-borne measure- 
ments that the 6300 A emission was approximately proportional to the maximum 
electron density in the F-region of the ionosphere. The dissociative recombina- 
tion theoiy was further advanced and developed by Lagos, Bellew, and Silverman 

[179] , Peterson, Van Zandt, and Norton [116], Holmes, Johnson, and Young 

[180] , Zipf [28], Wallace and McElroy [12], Ferguson, Fehsenfeld, Goldan, 
Schmeltekopf [80], Swider [27], and Okuda [181, 182]. 

Brandt [ 9] predicted that the 6300 A emission in the day glow would be 
approximately 50 kR, and Noxon and Goody [183] and Jarrett andHoey [132] 
obtained similar values. However, Zipf and Fastie [46] , Wallace and Nidey 
[184] , and Nagata, Tohmatsu, and Ogawa [185] got values of approximately 
5 kR. Noxon [10] reported observed values ranging from 5 to 50 kR with very 
large day-to-day variations. 

The first reported redline altitude intensity distribution from rocket- 
bdme measurements was by Gulledge and Packer [186] followed by Huruhata, 
Nakamura, and Tanabe [187] and Wallace and McElroy [12], The zenith 
intensity during the 1964 rocket flight reported on by Gulledge and Packer was 
approximately 22 R with no emission above 300 km with an intensity greater 
than 5 R. This result was consistent with the dissociative recombination theory 
in which the emitting layer moved upward during the night. The results of 
Huruhata, Nakamura, and Tanabe [187] were similar with an integrated inten- 
sity of less than 25 R and an emission layer from 200 to 320 km, peaking at an 
altitude between 270 and 280 km. Gulledge, Packer, Tilford, and Vanderslice 
[50] reported on the results of another rocket-borne measurement which showed 
that approximately 5 h after sunset the peak emission occurred at 246 km. 

The zenith intensity on the upleg was approximately 70 R, while the downleg 
intensity was approximately 90 R. The difference was attributed to an observed 
downward motion of the ionosphere, thereby lending more credence to the 
correlation between the redline emission and the ionosphere. 
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Eather [188] reported that the average night glow 6300 A emission was 
approximately 35 R north of 80° invariant latitude, then decreased to as low as 
5 R between 40 and 55° invariant latitude, and then increased to 30 R in the 15 
to 25° latitude region. These results were obtained from a tilting photometer 
flown on a Convair 990. The photometer system included a C 14 radioactive light 
source that had been calibrated at Fritz Peak Observatory by M. Gadsden. 

Schaeffer, Feldman, and Fastie [ 112] , based on the results of an analysis 
of a rocket flight at a solar zenith angle of 60°, reported that the measured 
6300 A emission above 100 km was approximately 1. 1 kR with 0. 29 kR due to 
photodissociation, 0.56 kR due to direct electron impact excitation, and 0.26 kR 
due to dissociative recombination. They did not claim high accuracy for their 
electron impact results but believed the photodissociation production rate was 
approximately equal to one-third of that expected from the solar flux values of 
Detwiler, Garrett, Purcell, andTousey [189]. However, their results were 
in good agreement with the solar flux values of Parkinson and Reeves [190] . 

Brasher and Hanson [ 191] have provided further evidence of the connec- 
tion between the ionosphere and the redline emission. They showed theoretically 
that a lifting of the F 2 peak by only 20 km reduced the peak night glow intensity 
by a factor of approximately three and in so doing confirmed the suggestion by 
Van Zandt and Peterson [ 192] that there is a correlation between the midnight 
enhancement of the 6300 A emission and the ionospheric drifts. 
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APPENDIX B 


LISTING OF COMPUTER PROGRAM OF THEORETICAL MODEL 
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7* 

T* 

n • 

1 u* 
11 * 
1’* 
17* 
14* 
15* 
IS* 
17* 

1 3* 
1 ?* 
?G* 
21* 
?.?• 
2 ’* 

2 4* 
2E* 
2S* 
27 * 

n* 

2** 

3U* 

SI* 

32* 

32* 

34* 

35* 

3 S • 
??• 
S3* 
39* 
4 o* 
41* 
42* 
43* 
44* 
45* 

4 fi* 
47* 
48* 
49* 

. 5U* 
51* 
5?* 

5 3* 
54* 
55* 
56* 
57* 


dimension ppnu 

DImE^cion 7FLUV<° 11 .COL U M 4 < 3.5®) #Z1 ▼** ) .E < 3*1 .X7( *> *TH{ ?) 

COMMON T AU (3) iT E^PO-» (41 • FLUX (45 > * S I? A«S (3. 4 5 )• AT *C1, * TMC , '». M -I'’*, 4 T 
1 M A f l TMS .ATMGAMC 3) .TONI ?G( 3) .E .FM TS €C C 50 L) .IN TE F S f f U 0 > .PEL 7 *F . J , 

2” IV F* C 4 A c 3. AL ATH.C *6N I. *L*»TH f C« 04, ** A OC. H* C 3 I • 01 C. 0?0. 03 P. 0* 

C PROGRAM FOR THF CA *- C IJL 4 T I? N DF F 30 LA »IRGtOW AT >®EfIFrre ALTITUDES. 

IN :5 
I 0" E 
IP - 7 

2UL FqpMA T(5«5. U .Fp. 1 .3E8. I .1 4 > 

201 c OPM AT ( 3E16-3) 

?L? F0PMATlElS.4flSXt5ElF.BJ 
2 l) 3 c 0 PM AT ( 1 3E 5 • 2 ) 

2L4 r CPMA T(SE 10. 2) 

2u5 FORM AT( 1 ? 1 5 J 

9G F0PKATC11H INPUT DATA/49H TIN F t r I ?(j ,C ON 1 ?D C 1 1 .CON 1 ?uf 21 .C CN 1 

l-’OtS), = /2( c lC.3flH* )«3(ElS«9*lHf)) 

99l F OP MA T ( 2 9H YE A R f DO .HR .MM . *HE T * r / l X . 4F 7 . 1 ,F * . 2) 

SUF FOPMAT(Fc.l. SF^.l) 

C NUM3ER DENSITIES TE MPO 7M> - TE HP 0 2( 4 1 


c DEFINITION OF CONST I*JFMS 

C I TEMPO ZC I) 

C 1 N 2 

C 2 02 

r 3 0 

r 4 oo 


C READ THE INPUT PARAMETERS 

00 17 Lrl .45 

17 PEAD( IN. ?u2J ^FLUX (L J. pTG«P c (3.L). S I c A»^ (2, LJ.SISAP^U f L) 

READ ( IN .5 U 4 ) 

READ (IN. 201) Fo N1 7j 4j m i C0N12Uf?). C0Nl?J(3) 

PE A 0 ( I N .20U) YEAR .T 120. ?4» X . ZM I N .PFL 7. Li AX 
10OU READ( IN. 5UFJT INF. DO.HPf MN 
PEA 3 f IN *?U3 ) ( E ( J ) . J= 1.39) 

WR ITE( 10. 90) T INF.T1 20. r 0 N1?U (1) . FONl 2u (2 ). CO N1 2 0 ( 3 J 
Tpl r 6. 28 T 18 5 35 
XCON =57.29577951 
EPS:23. 45 /XCON 

C FLONG: IS LONGITUHF OF HljNTSVlLLF 

r L0N3=-6S.F5 
PIT=1. 57074634 
PI=3. 14159768 

CALL ANSLE (HP. year. D". MN. X CO N. E PS . T P T # FLONG. T HFT A ) 

WRITE (10 .9 90) YEAR .DO. HR .4 4 .THETA 
00 14 1:1.5 
14 COLUMNCI . 1 ) =0. L 
85 CONTINUE 

C THF INITIAL VALUES 0 F THE ARRAYS - FLUXtt3.C0LuMNEI.JD. AND PRFVC 13 
C ARE NOW c STA9LlSHEn 

C THESE PARAMETERS ARE CALCULATED ONLY ON TF ANT THOSE WITH THF 

C PREFIX - ATM - ARE USED IN THF SU8R0'JTTNF H 0°AT M. AND APE SPEfIFIFP 
C IN COMMON 

n=TlNF-90U.D 

AR3:0/(75U. 0* 1.7 22 Z -4*0 ••?) 

ATMS=1 .5F-4*0.U291*EXP(-.5*AP G**?) 

A TMAr(TlNF-T12C)/TlNF 
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S3 C ATMC1 AND ATMT3 ARE IN CEMTIHFTERS* W HI L r A T HC2 I? IN *ILO*FTrftS 


59* 
5 G* 
151 * 
62* 
6 T* 

Sft* 

B5* 

So* 

67* 
5 5* 
S'** 
7G* 
71* 
72* 

7 Z* 
7ft* 
7 C * 
75* 
77* 
79* 
7«* 
90 * 
81* 
82* 
83 ' 
9ft* 

8 5* 
86 * 
B 7* 
89* 
B 9* 
90* 
91* 
92* 

9 3* 
9ft* 
95* 
96* 
97* 
99* 
99* 
QO* 
01 * 
02 * 
03* 
Oft* 
05* 
06* 
07* 
09* 
09( 
10 * 
11 * 
1?* 
13 * 


* TMC 1 = 1 20. GE 5 
*TMC?zSft76.77 
A TMC3=53 C S. 77E5 

ATNGAMf 1) =11 .1359*28. 0)/fATMS*TlNF> 

A T«SAHf2>=A TM3IM 11*1 32.0/29 .01 
AT*GAm( 3) = ATHG AN 12 ) * • 5 

J=t 

?A J) =?H AX 
72=ZIJ) 

T HET A =t het A/XCON 
CALL JACC 22! 

CAOC =9.UE- 1 1 
CA0A=5.QE-I1 
CA BN I =2 .UE- 1 2 
CAAFQ=ft.OE-ll 
ALBTH = 5.0E-U7 
A LA TH= t. OE -07 
00 5ft 45 1 = 1*3 

X7«,n = <ATHr2*77> /IN 2d l*l.E5l 

5Mft« TAUI I) =SQRT C .5*72 IT)* tC0S<THFTAI*CO; IT HET All! 

CALL ERFS(ERF) 

IFIPIT-T HEX A) 53 3* 5^3*535 
533 DC 599 1=1*3 

59 9 CH| I ) =| S QPT <P IT*X 2 II) *$ INI THET A) ) ) * F XP IT All II I **T AJ 1 1 1 ) *<1.*FRFII) ) 

50 TO 665 
535 00 60S 1=1*3 

6 06 CH(I)=IS0RT(PIT* XZ«I !*5lftlTHET«) M!»EXPITAL»| II*TAUCI)I*1 l. -ITRFII I I 
566 00 577 1=1.3 

5 77 COLUHNCl *J> = TEHPOZCI )*H2l I! *CH II) * 1 .E 5 
DO 66 L = 1 » 45 
FLUXIL!=0.0 
00 3333 1=1*3 

33 33 FLUXtL >= P LUXU>*SIPAPS( I *L ) * C) L UMN3 I * J ) 

FUilLIL! =7fLllX <L >*EXfM-FLUX CL) ) 

66 CONTINUE 

CALL CONSTI 

call emiss 

INTEGSI J) =FMIS5SCJ) 

wki tecio *iob)z?.i tempo zi r ) *i=i *ft> .mr sssui *ci J) .intfssi ji * 

1010.020.030. OP 
0 07 J=2 *500 

21 J) = 2<J-1) - 0EL2 

IF| 7IJl.LT. 2MlN)50 TO 21 
77 =2 1 J) 

CALL JACIZZ) 

DO 5ft 1=1*3 

XZII ) = |ATMC3* 22) /HZH >• l.E 5) 

5ft TAUI I) =SORT I .5*X7 II)* ICOSI THETA) )**2) 

CALL ERFSCERF) 

IFIPIT“THETA)53*53»55 
53 JO 59 1=1.3 

5 9 CHI I) =1 SORT CPIT*X2«I)*SINITHETA)) )*EXP IT Al) II)*T A'J II) ) *C1.*FRF|II ) 
SO TO 56 
55 00 50 1=1.3 

60 CHCI)=IS0RT(PIT*XZ1I)*SINITHETA)))!*EXPITAUII)* TAUIII9*! l.-FRFID) 
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14 

IF* 

IS* 

17* 

1?* 

I c * 

2 l* 

2 i* 
?*>• 

7 *• 
24* 

? t • 
*>S* 
•»?* 
73* 
23* 
70* 
71* 
37* 
*3* 

34 * 

3F* 
3f • 
77* 
33* 

4b* 

4 l* 
47* 

9 OF 


1* 

7* 

?» 

4* 

S* 


c * ^0 *7 in. 3 

5 7 COlUMMI *J) = TE*PO*M I 1*HZ< I1*M|T )• i .E 5 

*« •"’O $ L -1 * 45 
C LUX(L ):L. 

:70 3037 1=1.3 

3b ?* FLUXI L )= r Ll'X<L>* r t' fc * a S( I »L >*CDLUMN]| ! ,J> 

flhxil) = 7 flux «l >*exp i-fu»x iuj 
5 continue 
c*ll const I 
CALL ’Hies 
CALL INTEG 

UP I TEl IO.lOB) ZZ.( rrMPOZCT ) .1=1.4) .E«I rS5( J) .El J) .INT^SSI • 

1 0 1 0. 02 n . 0 3 0. OP 
7 CONTINUE 

UP IT El io. 103) 

30 TO ?2 

21 VI? IT cl 10. 104) 

2? ^CTO IGuO 

luU FORMAT! 1 x1H16HHEIGMt-*xFIU.T. IX7HKH.) 

1U T c CPM4 T( 1UX T 5HLCWER BOUNDARY 4 A * NOT B r FN REACHFQ) 

104 e O^MAT( IUX^SHMORHAl COMPLETION Of POH PUT A T 10 N) 

10F c ORf-ATI /Fx.ZMLr.IS.SXZHlr .12*? X2 4?= ,T 2 .* X 104^ A TA ERROR! 

1 b 3 c O'?MAT! IX. SH7Z = E7U . V IX 1 S HNU Me* R OEN? It T£ 5 / i x # ’ HN7 • lx El 1 .3. 2 X. *H. 
IP? . 1XE 1 1.3 . 7X.IH0 * 1XF1 1. 3 .2 X •2H0*' ,1 XE 1U 3/1 X.EHEMI5S .IXE 11. 3.2X » 
?IH*. IX El l .3, ?X. 6 41*1^5. tXFll .7X1 X . 3 HO 1 0, 1 X EU .3, ?X. 1 HO?P. 1 XEll .3* 
3?X*3HO?0 .IXE 1 l. 7 ,2 v . 2-1 OP . 1 XE 11. 3 XVI 
SU4 FO n H AT ( 7 R H 

l ) 

eno 

compilation: no orAPNOsnrs. 


eu?POl*TlNE JACI Z7) 

PEAL integs 
DIMENSION '(39) 

DIMENSION AL°HA IS) , El CGJ.OT! 51. Pf 7) f PTT (5 l.xM (J ) 

COMMON TAu< 3) .TEMPO ’‘I 4 ) .FL UX( 4 * 1 .? r G A 9 r I 3 .4 5 1 .A T MC 1 .AT mc 2 .A TM C3 .A T 
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7* 
9* 
9* 
10* 
II* 
1 ?• 
l»* 
1*1* 
I 5* 
13* 
17* 
1** 
19* 
2C* 
21* 
22 * 
2?* 
2*1* 
2*< 
?«• 
27* 
?9* 
2 ?* 
3b* 
31* 
3?* 
33* 
39* 
3f* 

3 5* 
37* 
39* 
39* 
4b* 
41* 

4 2* 

4 7* 
44* 
45* 
45* 
47* 
49* 
49* 
5b* 
51* 
52* 
53* 

5 4* 
5 C * 
55* 
57* 
59* 
59* 
50* 
51* 


X*A* ATMg, AT»6*M<3> * CON 1? 0 < 3) * E. FHISSE I'CD >* TNT W? !5GJ l.OELZ *F, J, 
2TINF *CAA*0 .ALATH .CABUT *AL9TH,CA0A *CA0C*M7(31 »D ID *0 2D *0 *D *C D 
T3b=193. 

7=7Z*l.UE-5 

5R A V =4.Su6S5E-3* Cl • 2/ AT NC3 ) 2 • 

ALPHA (41=0.0 
ALPHA! 7) =U.O 
A LPH4 ( ? ) = 0» 0 
ALPHA! 11 =U.U 
XHI 3lrl6* 

XM(2)=32.U 
T*M 11=25.0 

40LT?=3.7l432E-3/GP AV 
T=TlNF 

RH09Q = .34SE-9 
ALPHA (Sir-. 34 
ALPHA! S> ru.U 
A Vr6. 0?247 i 27 
EH 1) =?4.ul34 
£1(21=71.9954 
EX! 3) =1 5.9994 
El f 4 1 = 79. 948 
EH 5) =4.Cu?6 
El !6 Jr 1.00797 
*C 1) =24.92579 
3(?1 = -. 0740G66 
9( 3) =-.Oll 440 7 
9 (41=4.5 I1U3E-4 
9! 5) =-9.*i 495E-S 
9 I* 1 = 1. 0754 lE-« 

4( 7) =- F .97444 E- 7 
0N= . 7P 1 10 
90?=. 20955 
!?A = . 0093432 
9HE=S. 14715-6 
PK=S. 31432 

TX =371 .65 71 *.U5 14 805 *T-2 94 .353* *t*P (- .bU2JS?2*> *T 1 

A2=2.*CT-TXl/3.141 r 9265 

DIT! 6) =0. 

H=I 0 

EPSz.OUUl 

Tl=1.9*< TX-T9 0 1/35. 

T4=7.*!TX-T90 - ?!.* T1 *35 .*/ 4 .1 / < 3 5 .* *4 ) 

T3=-Tl/i 7.*35.**?)»4.*T4*35. /3. 

TZ=TX*T1* <7-125.1 3* l*-125. 1**3*T4* <7-125 .1**4 

IF .r^-lLU. I 4 3 *4 3 *40 
43 2 2 =2- 90 • 

E* = 3 ( II *Bl 2>*2?*9(7>*Z2**?*9 <41* -2* *7** ( C 1* T ?**4*9! F)***2** c 
1*9( 7) *72**6 

0=2 

7 (J CONTINUE 
Ar*G. 

FA=8( 1)*9<21 •< A-901 ♦* <31* < A-9Ul**2 *e<4l *<A-90 >••* *"<5 1* (A-»b )**4 
l*3m*(A-?0)**5 *8 < 7 1* ( A -901*45 
FA =FA*9. 9u 55 4/ ! < l • ♦ A/ 6. 756 76 6 F*7 1 • *2 1 

FArFA /( TX*T1*C A-125.1* T3* < A -12 5 . 1** 7 * T4* < A -12 4. 1 ** 41 
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ii i !■■■■■ ii iii i ii i n 1 1 in mi 


Mill! II II II nil i 


l 1 1 Bill I i I 


I I I I 


I INI I 



52* 



F0=9f t )♦ *17)* 10-90) *« 13)* tn-90)**2 *6 14) *10-90 )•** * 0 1S >*1P-9U>*j*4 


si* 



l*91F)*f0-9QJ**5 *9171*10-901**6 


-'4* 



FD=FD*4. Ru56*>/1 <1 . *0/ 6.356T6G E*7 >•♦* > 


ss* 



r DzFD/CTX*Tl*10-125.)* 13*10-125.) *4 3 * T 4 * 1 0 - 12 E. > * * 4 ) 


S3* 

r. 


SR34. SIMPSONS RUL E QUADRATURE - 3.F.KUNCIR 

SR34UL 

S 7 * 



definitions - 

SR34W 

S9* 

r 


* z LOWER L I w IT OF INTF3PATI0N 

SRQ4CI 

S3* 

r 


0 r UPPER LIMIT D- I N TFGR A TI ON 


7u* 

c 


FUNC = INTEGRANO FUNCT ION SU^PR OGR AH 

SRQ4Ui 

71* 

c 


EPS z RELATIVE EPPOP CONVERGENCE T9 I T r R ID N 

5 R P 401 

72* 

r 


0 r RESULT OF INTEGRATION 

SRC40( 

77* 

C 


N r NUH3ER OF INTEGRATIONS RE Q Ul p EH TD cjnq p. 

5R04QI 

7 *4 • 




69Q4UI 

7** 



N = 0 

SR04Ui 

7'* 



RR Ev -0 . 

SRQ4CI 

77* 



SONEz (0-A )*(FA* r D) /2. 


7 ? * 

71 


N zN* 1 


79* 



TF 1N-M) 7 2 » 72 »7 5 


*a* 

V 


** I N T Z2**N 


91* 



S TUOZO. 

SRQ4L1 

9 2* 



OELzlD-Al/PLOAT 1NINT) 


97* 



00 73 Izl.NlNT,? 


99* 



X=A*OEL* C LOAT IT) 

SRQ4GC 

9T* 



Fx=5( 1 1*3 1 21*1 X-®0)*»( 31*1 X-RL1**2*RI 41*1 X-901 *• T + E 1 5 ) • 1 X- 4U) * • 4 


93* 



1*91 5)*< X-90) *** *9 IT) * lx-90)**S 


9 7* 



r X=FX*9. PUE65 n < !.♦ X/S. 3S6766E* 3)**?) 


9°* 



FX zFt/| TX *T1* (X-12^.1 +T3* (X-l 2S.)**3 *T4 *fX- 1 2 5 .) **4 ) 


PB* 

75 


c TWOz S TWO* r X 


4C* 



CUR =S0NE* 4 .* OEL * c TWO 

SR 04 Ul 

* 1 * 



IF lEP<“*ft 3 Si CUP ) -A“ SI"* UP-PPEVH 7 4', 75 ,7? 


9->* 

74 


PO r V zCllR 


9 7* 



SCNEr c cCNE*C UR 1 /4. 

SRG U(j| 

94* 



GO TO 71 


9 C * 

7* 


PrClIR/* 


9S* 



IF 1 Z- luO.) 44* 75, 44 


97. 

4* 

IF {0-1 LU. 1 76 *55,76 


99* 


7* 

PENS zR HO 9U*T 9 u*FM*FXP 1-9 f FK ) / <T? * 2 4 . 9? S7 9 1 


90 * 



0 L r XL0310ENS) 


Ul* 



oApr AV*OF,n^/FM 


Ul* 



A Nr XL0 31 3 N * Z M * P t P /?B . 3 6 1 


u?* 



IAZXLOG 13 A*E«*P *R/?9. 9S) 


\jXm 



AHEzXL0G(3HE*EM*PAP/29.?SI 


U4* 



A0=XL0G(7.*PAR*11.-Fw/23.95)) 


05* 



AO?zXLOGlPAP*C r M*l 1 . + 0 0 ?> /24 . 9 S - 1 . 1 ) 


U = * 



AHz-O. 


U 7 * 



TE M PO ?C llzAN 


Q 9* 



TEMP0Z1 2) =*0? 


U9* 



TEMP07<3)rA0 


1 U* 



003331 1=1.3 


11 * 

7371 

H ?1 1 ) = (30L t 2*T7I / XMl I) 


1 *>' 



RETURN 


17 


40 

73- I Ob. 


14* 



TZTzIX*T l* 1Z 3-125.) *T3* (7 3-l*?s)**3*f 4*17 3- 12 5 )-**4 


IT* 



7H3z3ll)*6( 2)*TU.* 9 C 7) * I DU.** 1 4 ) * I Ouu. ♦» 1 E) • IL. * • 4 •* E 1 6) • 1L •» *f • 


1 -• 



1 *91 7) * 1L .**5. 


17* 



Dz 1 00 . 
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1? 

1°* 

2 u* 

21* 

2 5 * 
74* 
25 * 
23 * 
27* 
? 9 * ; 
??* 
30* 
31* 
3?* 

IT* 

34# 

T r * 

35 * 

37* 

33* 

39* 

40* 

41* 

4?* 

45' 

44* 

45* 

4** 

47* , 

4 9* 

4 2* 
SC* 
M* 
5?* 
s J# 

54* 

55* 

S3* 

57* 

S 3 * 

* ?* 
6 u* 
61* 
6-»* 
5 ’* 
54 * 
c5* 
56* 
57* 
59* 
59* 
7U* 
71* 
7?* 
77* 


50 TO 70 

55 ?ENl=PH09U *T9Ci * Z* 3*E X a ( -2 /FK ) /( T*>7* 25U P2*79 » 

P4P = *V*0 f :n1/7M? 

OX ( 1 ) =QN* ZS 3*PAP/?«. ®6 

p If .?) =*A9* C2M3* (1 •♦002)/ 7 8- 95-1.) 

QI<3>=?»*Pt**< t.-Z*2/29.«>5) 

0I< 4) r9**Z^3*P4R/2*«96 
DI(5)=Grt5*7H3*PAP/??.96 
r Ff7-l’5.) 56* c 5. 90 
**6 CONTINUE 
M=lUO. 

F4l = *»9U'S5/CI 1* * A 1/6. *5 €76 6F ♦*)**?) 

FA1=FA1/<TX*T1* <A1-125.)*T3* t A1-12*.)**3*T4.*1A1-1 25.)**4) 

01=7 

F01=9.4oS65/ f f 1 .*01/6.355 766 E Ml* *21 
IF (01-125. ) 45 *45 #50 

4 5 FClrFDl/ ITX»T1* f 01 - 17 5 . } ♦ T 3* 101-1?*.) • *3 4 • < 01 -J 25 . 1 • *4 ) 

SO TO 51 

Su F01=FR1/<TX*A2* AT AM It 1* (01-125.1* <1. *4 .5 E- 6 MCI -t 25.) **2.S I/A2 )) 
T2rTX*A2*ATAN( TIM 125. 1 * « I .♦ 4 . 5E -F* C 7 -1 25 • ) • *2. 5 1 / A2) 

N = 0 
P©5 V=U 

SON£=( Ol-AllMFAl ♦P01)/2. 

91 NtM* l 

IF f N-*1) 9’* 42# 9S 
*2 NlNT=?**N 

STWOrC. 

DEL- C 0 l *4 1 1 / c LOfi T(WINT) 

00 93 I = 1 #N INT # 2 
Xl = AlOEL* ff LOAT<I ) 

FX J=9.4LS6S/ f 11 .*X 1/6.755756 E*3>**?> 

I F f XI -1 25 • ) 46 #46 #52 

4 5 FX 1 rFXl/ (TX*T1* IXl-17 5.) *T3* 1X1-1 2 C .) * *3 ♦ T 4*fXl -1 25.) **4 ) 

30 TO 93 

r 2 FX1=FXI/(TX»42*AT«N<T1* (X 1-1 25 . ) * <1 . +4 .5 E- 6 MX 1 -1 25 . I **2 .S 1 / A2 )) 

93 5TtfO=STWO*F XI 

CUP=S0NE*4.* OEL *STWO 

IF (EP6*AB 5CCU») -M 5CT UR -PPEV) ) 4 4 #95 #9 5 
34 PR EV =CUR 

*0NEr<cQNEMUP)/4. 

50 TO 41 
45 P = C UR / 3* 

00 41 1=1# 5 

DIr(I)r0Kl 1*1 r?3/rz>**f I.*»L»H*m }*FyPt -FI II )*R/ fk ) 

4 1 CONTINUE 

T E*PO ZC 1 ) =3 I TC 1) 

TE«P0Z<2)=0IT(2) 

TEMPOZf 3)=3ITC 3) 

OEMS =0 
00 42 1= 1 #6 

OENS=OENS*EHI)*niTfn/Av 
42 CONTINUE 

EMrOENS* Ay/ <OlT U 1 +OTT (2) *011 IT) ♦ OTT <4) ♦0IH5 )*OTT C5 ) ) 

?L= XLOSOENSI 
AN rXLOG (Oil Cl ) ) 

*o2=xlosoi tern 
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74 

75* 

7 5* 

77* 

73* 

7° * 

3b* 

31* 

3?* 

3** 

33* 

35* 

35* 

37* 

*3* 

3 3* 

7C* 

31* 

4** 

9?* 
34* 
95* 
43* 
9 7* 
99* 
99* 

0J% 

<Jl* 

0?* 

b T * 

U4* 

ti c * 

05* 

07* 

0 9* 
J9* 
lu* 
11* 
1 ?* 
1?* 
19* 
15* 
15* 
17* 
19* 

1 4* 
2U* 
21* 
2?* 
27* 
24* 

2 5* 


AO =XLOGCOlT C3> ) 

• A r XL03C D I T( 4 ) ) 

AHt!zXLOG«DlT (S> > 

I c t Z-500. 1 47 *4A .44 

47 OITC 6) =1C.** <-6> 

49 Hr XLC3CDITC6)) 

an =A*AX1C-U** AN) 

A C2r A HA XI ( -U. *10?) 

AO =ANAxU-0.* AO) 

I A :I9M1(-U. . A A ) 

AHE = AMAX 1 <-U.* AHE) 

AH r I XIC -0. * AH) 

DO 3*32 1 = 1# 7 

.337? HZlI ) = (90LTZ*T7) /XMC T> 

° ETUPN 

9U r =rxa?*A TANC Tl* 375. •< l. ♦ 4 . 5 F- F* 3 75 .** *> . 5 ) ft 2) 

PI( 6) =lb.**C ?3.13-’'4.4*XL0<MS) *5.5*X LOGtS) *xLOGC e ' ) ) 

A IrBUD. 

IFCZ-400.) 44, CQ, 60 
44 A 1 = 2 

50 FAl=9.9uSSF/ Ctl • ♦ Al /S . 3 53766 F *7 )• ♦ 2) 

FA!rPAl/< TX*A2*ATAN( Tl» U 1-125 .)• C 1 •« 4« $*■ -6* U l- 125. I* *2. M/A?) ) 

01 =z 

I p ( Z-500. ) 6l*F2*S? 

FI 01=500. 

6? FDl=9.B0665/(C 1.+D1/6* 3567S6E*7)**?) 

F01=F01/<TX*A2* AT AN CT 1* CPI- 125*)* Cl, *4 .5 €-6 *101 -12*?.) **2.5 >/A2 )) 
N=U 

PR EV =0 

50NE= (? l-A 1 ) * ( c A l * r D 1) /2m 
31 N = N* 1 

IF ( n-m) 9 ? *3 2*95 

**2 NlNT=2**N 

r TWO = l>. 

OEL=( 01- AD/FLO AT <NINT) 

DO 97 1 = 1 *N IN T *2 
XI =A1*0EL*FL0AT II) 

c Xl = 9. 9 0655 /<< l.+Xl/F. 35 6 76 6E* 3 > * * 2 > 

FX 1=FX 1/ CTX*A2*AT AN IT 1* (X 1- 1 2 5 • ) • C 1 - *4 .5 5 * <Xl -1 2 5 . ) * *2 .9 ) / A2 ) ) 

0 3 S TWO= S T WO* c X t 

CUP =5 ONE* 4 •* DEL **7 WO 

IF (EP«**A3 5CCUP) -•* SIT UR -PR E V) ) *4,45 .9 5 
44 PPEV=CUR 

F0NEr(C0VE*CuP>/4. 

30 TO 91 
»5 o = rUR/7. 

▼7=TX*A2* ATANITl* 1^-125.)* Cl. *4.5E-K* C7-125,)**? .5 )/A2 ) 
IFCZ-5UU.) 6 3.64.64 
5 7 P=-R 

64 DIT(6)=0I C5 )*CF/TZ)*FXPC -El C6)*7 /FK ) 

00 TO 56 
-END 


i o- compilation: 


NO DIAGNOSTICS. 


105 



i* 

2* 

?• 

4* 

E • 
7* 
• • 
1* 

IL» 
11 * 
12 * 


SUBROUTINE INTF6 
REAL I N Tt 3 S 
OIhENSION 

COK*ON TAUC 31 »T r MPO *M H) t FL UXl 4 4) *ST FA? S « 3 1 4 « I • * TM- 1 ,AT NC 2 ,A TH C 3 ,A T 
lH A v ATRS 9 ATMGiHl3) # TON t?Uf EtEHISS? <*U2 ) f INTERS < 50J ) • DELZ. F, J t 

2TINF »C AA ? 0 •ALATH *f ABNI 9 AL§ TH t CAO A t CAOr f H7«3l *0 10 *0 20 ,0 3? *0® 

J J - J * I 

• EMlSS = <tPl SSS<JJl*EPlSSS(Jn*.S 

C0L£Hl=AFHlSS*PEL7 

INTESSI J) = lNTE3sOJl*C0LF4r 

RETURN 

END 


1* SU°ROUT INE ANGLC<HR,YR9 nAYNO.NIN.xCON, EPS.T°I. FLONR, THETA 1 

2* PE 1 L jDtMIN 

3* JD=24398SS.*(YR-1 46 9 •)* 3*5 •♦BAY NO 

4* 170 TEmP=| JO-24 1SL2LO/3FF??. 

5* TGnr9«*S4U9433*7E" p * < 36 00U . 76 854 *3 .4 7U RE-04 *T EMP I 

F* 171 I F f T30-3SU. Il73.l73fl72 

?• 17? TGOrTGO-160. 

5* SO TD 171 

7* 173 GNT=HR* 60. IN 

U* 17F TH3 = T30*3HT*.2Fl;FB447 

t 2* 174 IF<THG-36J.) 1 74t 1 73tl 77 

12* 177 THS= TNS-75U. 

I** GO TO 17«5 

14* 178 AN3rTH3* r L0NS 

1«* 130 IF( ANG-3SU013U*l«a*l 35 

ie* 1R5 ANSrAN3-36U. 

17* GC TO 180 

14* 190 ANSrANG/XCON 

1 <?* FL *HS =«0l 7?u 3* C Jn.r4 3 5 33 a .) *.U33 ,; *S INI .Jl 7 20 3 * 

20* 11 JO -?4 *5* 3° •>)-!• 41 

’ 1* IF( FL A“S 11 35 *191*1 91 

>2* 191 IFIFLANS-TPI) 1*9 »1*3 .192 

’?• 13? Fl AMS=FLAHS-TPI 

>4* 30 TO 191 

?3* 1 9S IF(FLANS*TPl)196f l*9t 199 

?S* 196 F|_ANSr p LAMS* TP I 

27* GO TO 196 

29* 1 98 CONTINUE 

21* DL SUN r AS IN(S INIEPS)** IN (FLAMS)) 

30* ALSUNrASlNl Tl N ( D LSUN 1 / TAN C E * SI J 

31* C PUT ALSUN AND FLA*S IN SANE QUAfV 

32* AtSUNrA 3 SI « LSUNJ* xroN 

3 3* TEMP=A9S (FL ANSI *X CON 

34* IF | UNO -*CU l?3(J»?30*?ULi 

35* 200 IFITEMP-I 3D.)205f 205t 210 

J** ?CE ALSUNrlRG. -ALSUN 

37* GO TO 230 

35* ?!0 IF|TSHP-?7C.)215t2l5t22b 

39* 21 S ALSUN = ALSUN*1B0. 

40* 30 TO 230 

419 220 ALSUN =360 •- ALSUN 

42* ? 30 IF(FLA4S1235»240 924L 

43* 23S ALSUNr-ALSUN 

4 4* 240 A LSlfNr A LSUN /XCON 

4 f ♦ HSi|N=ANG- ALSUN 

46* HRirHSUN 

47* OS=OLSUN 

• 8* C THE TA rCOsI 34 . 6 / XCON J*COS( n 5 1* CO SI M.9 A >♦ <1K( 34. XCCN>* SIMPS ) 

43* T HET A r ACOS CCT HFT A) 

SO* THE TA r THE TA *XCON 

SI* RETURN 

2* END 
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u function xloscxi 

?• XLOGrALOGlOlXI 

3* PE TURN 

4* END 

d c compilation: no diagnostics. 


lft SUBROUTINE CONSTl 

2ft REAL INTEGS 

T* DIMENSION 5139) 

4* COMMON T AH (31 • TEMPO? < 4) • FLUX 145) * S IG APS 13? 45 I* AT MCI* ATMC?* AtMCX. AT 

c* IMA •! TKF tATMGAMl 3) *rONl20( 3) *E*E9 Tss*f 5UCi> *IN Tt G S < 5UL> *DSL Z*F*J* 

8* 2T INF* CA AFO* At AT H* C A8N I* AL*TH* CAOA. CAOr, H* 131* 01 0* 0? 0* 03 C* 0® 

7* OrECJl 

9ft RP zQ • 30 

9ft PRr 1 5 

10* PRR sO • 1 5 

lift p=U0 

12ft SUMlrO.U 

13ft DO 51 L- 1*4 5 

14ft DISSOCrSlGAes <2,L >*FLUX 1L) *T EMPO? 17) 

!*• 51 f UM l: SUM 1* D 1 5S0C 

15* 0 1 OzRftSUM 1 

17* 0P = 1 • P5*Q * *? • ) / ( 0* T r 4P0?l 1 )*CA 3N I /A19 th* TE HP 0 Z( ?1 *CAA* 3/AL AT Hi 

19* 0?0=RR*CAAF0*TFMP07 12 I *0P 

14* 03D=RRPftCARNl*TEMP07C 1)*0P 

?0* TEMP02I 4) rOl 0*0? n 

>1* RETURN 

22* ENO 

) c compilation: no diagnostic*. 


1* SUBROUTINE EMISS 

’* PEAL INTEGS 

3* DIMENSION El 39) 

4 * COMMON T AU f3)*TEMP0?(4) *PLUX 1*51 • 5 IF APS 1 3* 44 )« AT MCI* A TMC?* ATNC3, AT 

5* IMA ,4 TMS *A TM6AM I 3) *rONl20t 31 .E.FUSSSl 5U0I »XNTE5S(50U) *OELZ*F*J* 

5ft 2 T lNF # C4AFo* ALATH, CA3NI. ALPTH*CAOi*CA OC* H? I 3 )* 01 0* 0?0* 03 P* OP 

7* » = 5.9E-9*TEMP0?«4> 

Bft FMISSSlJ):A/l9«l F -3*CA0r*TEMP07fl)ftC AO AftTEMiPO? 121 I 

9* RETURN 

10ft END 

3 OF COMPILATION! NO DIAGNOSTICS. 
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I II II 


1* SUBROUTINE ERF^IERF) 

2* DIMENSION EPP(5> 

7* COHN ON T A J 13) 

4* * lr.?F4829S9 

S* A2=.?34496?4 

6* *5ri. 4214137 

7* *4r-U4S11S? 

8* »Fr 1.LS14UF4 

3 * *=. 3274911 

1I> 00 1 1 = 1,3 

u* os=i./< u^*T*u(m 

12# 1 FR r (I|rl»-( tC(A5*Q**A4)*G5*A3)*9c.,i ?I*0 C *A1 lO r' r XP|TA'Jf T 1 *TA Ut I ) ) 

17* RETURN 

14* EN3 

3 o c compilation: no OH~N0STlCS. 
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APPENDIX C 


OBSERVATIONAL DATA 


The airglow equipment was operated on approximately 15 nights during 
the period November 10, 1971, through September 1, 1972; however, equipment 
status, weather conditions, and the requirement for a moonless condition 
limited to 5 the number of nights on which observations were made simultaneously 
by the C-4 ionosonde, Fabry-Perot interferometer, and the turret photometer. 
The temperature histories of the nights of January 7, 1972, and February 7, 1972, 
both contain very low temperatures that may be due to contamination of the 
Fabry-Perot 6300 A OI emission observations by OH emissions. This is dis- 
cussed in more detail in Appendix D, Electron density-height profiles, required 
for a detailed analysis, were not available on the nights of November 17, 1971, 
and February 4, 1972, due to equipment malfunctions. The night of January 6, 
1972, was selected for analysis because the records were somewhat more 
complete on that night; however, the data from the remaining four nights are 
provided here for comparison. The intensity time histories of all five are quite 
similar; however, the temperature histories of all five show large variations in 
short time periods with extremely low temperatures being observed near 
midnight on January 7, 1972. These temperatures are comparable to tempera- 
tures derived on other nights when observations were begun long after foe sunset 
period; however, to date no explanation has been found. 

Although the peak electron number density-time histories are similar, 
the actual number densities on January 6, 1972, are smaller foan those 
observed on the other four nights. 

It is believed that these data will support the contention that foe data of 
January 6, 1972, are not atypical of the airgLow measurements over Huntsville, 
Alabama, and, therefore, the conclusions reached on the basis of the analysis 
of fois single night are truly representative of foe majority of atmospheric 
conditions and not merely fortuitous. 

Figures C-l through C-7 present the observational data. 




Figure C-l. Observational data for November 17, 1971. 
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Figure C -2. Observational data for J anuary 7, 1972. 
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APPENDIX D 


CONTAMINATION OF THE 6300 A 01 EMISSION BY 
HYDROXYL EMI SSI ONS 


The dominant emissions in the night sky come from radiative transitions 
of the OH radical. The intensity and temperature of the OH( 9-3) band, which 
has its origin at 6356 A [ 16] , is important to an analysis of the 6300 A OI 
emission because the OI line is embedded in this band system. Therefore, 
since both the turret photometer and Fabry-Perot interferometer filters used 
to obtain observational data in this analysis are of finite thickness, they trans- 
mit some of the OH (9-3) band emission along with the 6300 A OI emission. 

The problem is to determine the magnitude of this contaminating emission. 


A. Photochemistry 

The detailed analysis of the OH band system by Meinel [193] showed 
that there was an absence of OH band emissions originating from vibrational 
levels with v > 10. Bates and Nicolet [194] and Herzberg [195] suggested 
that the source of the observed radiation might be hydroxyl radicals formed 
by the following mechanism, since the excitation energy of the v = 9 level in the 
electronic grdund state of OH is 3. 32 eV: 


0 3 + H — OH + O z + 3. 32 eV 


(D-l) 


Hydrogen was proposed to be regenerated by 


OH + O — H + 0 2 


(D-2) 


The spectral characteristics of the OH night glow have been reproduced in the 
laboratory by the ozone mechanism, equation (D-l) [196] . A photochemical 
model by Bates and Nicolet [194] indicated that this reaction, equation (D-l) , 
was capable of producing the observed emission rate. 
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Krassovsky and Lukash enia [197J suggested the following alternative 
mechanism involving vibrationally excited oxygen molecules: 

0 2 * (v > 27) + H -* OH*(v= 9) + 0 , (D-3) 


with the hydrogen being regenerated by equation (D-2) . Excitation of OH to the 
v > 10 vibrational levels, which required 0 2 * ( v > 26) , was precluded by the 
following: 


0 2 * (v > 26) + 0 2 — *■ O + 0 3 


(D-4) 


This mechanism was disputed by Bates and Moiseiwitsch [198, 199]. 

Although no conclusive evidence supporting either of these mechanisms, 
equations (D-l) or (D-3) , is in the literature, the general consensus seems to 
support the ozone mechanism, equation (D-l) . 

Table D-l contains a more detailed list of the reactions, with their rate 
coefficients, that might be involved in the OH chemistry. The principal reac- 
tions involved in the OH chemistry are shown in Figure D-l with reactions of 
equations (D-l) and (D-2) currently believed to be the most important. 

Table D-2 lists the reactions, with their rate coefficients, currently 
believed to be involved in the atomic hydrogen chemistry, while Figure D-2 
shows the principal reactions. Table D-3 and Figure D-3 contain the same 
data on the ozone chemistry. These tables and figures are from Reference 34. 

1. THEORETICAL PHOTOCHEMICAL MODELS 

Temporal variations of the OH intensity are of great interest since they 
are closely related to the chemistry of the D-region of the ionosphere. Bates 
and Nicolet [194] were the first to attempt theoretical calculations of the night- 
time kinetics in an oxygen hydrogen atmosphere. Hunt [104, 122] , Gattinger 
[200] , and Hesstvedt [201] have also theoretically investigated diurnal varia- 
tions in the OH radical production and loss rates and their height profiles. 
Hunt’s [122] results show that the OH production rate maximum occurs at 
approximately 40 km altitude at noon and at approximately 80 km at midnight 
and while the total columnar production rate decreases by approximately 
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TABLE D-l. CHEMICAL KINETICS OF OH 


Reactions Forming 

Reactions Removing 

O' + H^-OH + e 1 X 10‘ 10±2 

' 

oh + h + m-^h 2 o+m 1 x 10‘ 3I±1 

HjO + hi>-K)H + H X < 2390A 


H 2 0 2 + hv -*> HO + HO X < 5650A 

OH + O -+H + 0 2 (5±2) X to* 1 1 

O^D) + H 2 0 -* OH + OH 1 x 10* 11 

OH + H->0 + H 2 9x 10‘ 12±0 - 3 e * 380 °/ T 

H + 0 2 ->OH + 0 4 x io' 10±0 - 3 e ' 840 °/ T 

OH + H 2 -►HjO + H 6.3 x 10‘ 11±0,3 e - 27S °/ T 

H + 0 3 -► OH + 0 2 1 x 10" 1 1±2 e ' 500/T 

OH + OH -^0 + 0 1.4 x 10' 11±0 - 3 e ' S00/T 

o + ho 2 -+oh + o 2 1 x 10' 11 

OH + Oj ->02 + H0 2 1 .5 x 10" 1 1 T 0-5 e ' 4 ^ RT 

H + H0 2 -*0H + OH 2 x 10‘ 10±1 e ' 500/T 

OH + H0 2 -► H 2 0 + O 2 2x10' U±1 

0 + H 2 0 2 -►OH + H0 2 1 x 10' 12±2 e - 200 °/ T 

OH + H 2 0 2 -►H 2 0 + H0 2 3 x 10' U±1 e * 500/T 

H0 2 + 0 3 -+OH + 0 2 + 0 2 1 x 10‘ 14±2 

OH + OH + M -►H 2 0 2 + M 3.3 x 10‘ 30±0 - 5 (T/300)' 1 

0 + H + M -►OH + M 1 x 10' 33±1 

1 





Figure D-l. Chemical kinetics of OH principal reactions. 


60 percent between noon and midnight, the integrated intensity remains approxi- 
mately the same. This is because the higher production rate in the daytime 
occurs at a lower altitude where the excited state is collisionally deactivated 
before it can radiate. Intensity-height profile results from Hunt's [122] 
theoretical calculations also agree very well with the rocket observations of 
Packer [202], 

OH intensity temporal variations in the night glow can only be due to 
changes in the concentration of hydrogen or ozone. These changes may be due 
either to changes in the altitude of the emitting layer or to changes in the 
reaction rate coefficients due to changes in the ambient temperature. According 
to Sivjee [ 144] , theoretical calculations of ozone and hydrogen concentrations 
at a fixed altitude based on chemical equilibrium concepts show long rise and 
fall times, on the order of 6 to 12 h. These results are in serious conflict with 
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TABLE D-2. CHEMICAL KINETICS OF H 


Reactions Forming 


Reactions Removing 


H 2 0 + hv -*■ H + OH X < 2390 A 


H + Q'-*OH + e 1 x 10" 10±2 


O + OH -* H + 0 2 5±2x 10‘ : 


H + 0 2 + M -*■ H0 2 + M 5±2 x 10" 32 (T/300)" 2<3± ^' 3 


0H + H 2 ->H + H 2 0 6.3 x io' 11±0 - 3 e ' 275 °/ T 


H + 0 + M-*OH + M 1 x 10" 33±1 


H + H + M -*H 2 + M 1 x io ' 32±0 - 5 (T/300)' 1 


H + OH + M ->H 2 0 + M 1 x 10" 31±1 (T/300)" 2 


H + 0H->H 2 + 09x 10' 12±0 ' 3 e ' 380 °/ T 


H + 0 3 -> OH + 0 2 1 x 10" 1 1±2 e " 500/T 


H + 0 3 -+0 + H0 2 1 x 10" 12±2 e 12 °/ T 


H + H0 2 ->H 2 0+01 Xl0' 12±2 e" 500/T 


H + HO, -*• OH + OH 2 x 10' 10±1 e " 500 / T 


H + H0 2 ->H 2 + 0 2 6 x 10" U±1 e " 500/T 


H + H 2 0 2 ->-H 2 + H0 2 1 x 10" 12±2 e - 180 °/ T 


H + O, -+0 + 0H4x 10 -10±0.3 e -4800/T 


2 : 






Figure D-2. Chemical kinetics of H principal reactions. 


observational data that show short-lived variations [144] , on the order of 
minutes to hours; therefore, constituent number density variations at a fixed 
altitude do not appear to be the cause of the temporal variations. 

Sivjee [144] concludes that the short-lived variations in OH intensity 
are due to changes in the altitude of the emitting layer with the higher intensities 
occurring when the layer moves up to the altitude of the atomic oxygen number 
density peak near 95 to 97 km. These higher intensities are primarily a com- 
bined result of an increase in the atomic oxygen number density that controls 
the Amount of 0 3 available to react with atomic hydrogen through 


O + 0 3 + M - *■ 0 3 + M 


(D-5) 
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TABLE D-3. CHEMICAL KINETICS OF O a 


Reactions Forming 

Reactions Removing 

°3 + N 2 + _>0 3 + N 2 2±1 - 5 x 10 ' 7 (T/300)‘ 1±0 - 5 

0 3 + e -> 0 3 ' + hn 1 x 10' 17±2 

0 3 ‘ + 0 + ^ O 3 + O 2±1.5 x 10 ' 7 (T/300)' 110 - 5 

0 3 + e -+0* + 0 2 3 x t0 * 12±1 - 5 

0 3 ‘ + NO + -»-0 3 + NO 2±1.5 x 10 ' 7 (T/300)' 1105 

°3 + 0 2 ‘ -»0 3 ' + O, 4.0 ±1.3 x 10 ' 10 
0 3 + 0'->0 3 '+0 5.3±2x 10 ' 10 

0 3 " + hi>->- 0 3 + e 6 ± 2 x 10' 2 

°3 + °2 _> ‘ 0 + 0 2 + Oj 1.7 x 10 ' 12±0-3 e ' 143 °/ T 

0 2 ' + 0 ->-0 3 + e 2.5±1.5 x 10 ‘ 10 

o 3 + o 2 ^S) ->-o 2 + o 2 6 x 10 13±1 

O' + Oj ( J A) -*• 0 3 + e 5 x 10' 15 

Oj + O -^0 2 + 0 2 1.4 ±0.3 x 10 12 e (' 1 500±200)/T 

O' + 0 2 ( 1 S)^-O 3 + e5xl0 ' 15 

0 3 + N ->NO + 0 2 5.7 x io ' 13±0 - 3 

0 3 ' + N0 2 ^-0 3 + N0 2 7 x 10 ' 10±2 

Oj + NO -* 0 2 + N0 2 9.5 ± 1 x 10 ' 13 e (* 1 300±100)/T 
0 3 + N0 2 ->-0 2 + N0 3 2.6 x 10 ' 13±5 e (-7000±600)/T 
0 3 + H ->-OH + 0 2 1 x 10 ' 11±2 e ' 50 °/ T 

° 2 + + N0 2 "*■ °3 + N0+ 1x10-1 1±2 

0 3 + H ->0 + H0 2 1 x 10 ' 12±2 e - 12 °/ T 

0+0 2 + M->-0 3 + M 5.5±2 x 10 " 34 (T/300)' 2 ' 5±03 

0 3 + 0H-*-H0 2 + 0 2 1.5 x 10' n T°' 5 e ' 4|/RT 
0 3 + H0 2 ->OH + 0 2 + 0 2 1 x 10 ' 14±2 
0 3 + 0( 1 D) ->-0 2 + 0 2 7 x 10 ' 12±1 

°2 + °2 + °3 1 x lo1 1±2 

0 3 + NO + -> N0 2 + + 0 2 1 x 10' 13±2 


and an increase in the reaction rate coefficient for equation (D-l) due to an 
increase in the ambient temperature, although the atomic hydrogen number 
density does increase somewhat with height at these altitudes. His calculations 
are in excellent agreement with the rocket observations of Packer [ 202] and 
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+ M -*■ 0 3 


+ M 



Figure D~3. Chemical kinetics of 0 3 principal reactions. 


Baker and Waddoups [ 203] . When the emitting layer moves upward, it moves 
where the ambient temperature is higher; thus, his calculations support 
the observational data of Krassovsky et al. [204] andHuruhata [205] that there 
is a dependence of OH intensity on rotational temperature. 

If Sivjee's [144] conclusions are correct, a dynamical atmosphere at 
D-region altitudes, approximately 80 to 100 km, might provide a possible 
explanation for the extremely low and rapidly changing temperatures derived 
from Fabry-Perot interferometer measurements (Fig. C-l through C-5) of 
the 6300 A OI emission. The configuration of the F-P interferometer used in 
this study is such that any OH (9-3) band contamination is added to the 6300 A 
OI emission line profile at the peak. The addition of any OH emission that is 
at a lower temperature (approximately 150 to 250 K) causes a reduction in the 
derived temperature with the reduction being almost linearly related to the 
amount of OH intensity present as well as the rotational temperature of the 
emitting OH radicals. For example, if the observed emission consists of 
80 percent 6300 A photons from *D atoms at a temperature of 1000 K and 20 
percent 6306. 837 A photons from OH radicals at a temperature of 200 K, the 
data reduction technique used herein would derive a temperature of approxi- 
mately 800 K for the emission. Additional observational data are required to 
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separate the possible two components that may have resulted in the low tempera- 
ture data shown in Figures C-l through C-5; however, it appears from Sivjee’s 
[144] work that vertical motions on the order of 10 to 15 km due to vertical 
circulations or acoustic or internal gravity wave propagation could very well 
explain the observed abrupt, short-lived changes in intensity and/or temperature. 


B. Observations of Intensity and Rotational Temperature 

The original work on the OH emissions was done by Meinel [ 193] who 
found that for the time period from about 1. 5 h after sunset to approximately 
1. 5 h before sunrise, approximately 50 percent of the 6300 A OI emission 
measured with a broad bandpass filter was really due to the OH (9-3) band 
system. Chamberlain and Smith [205aJ , on the basis of theoretical calculations, 
predicted that the intensity of the OH(9-3) band should be approximately 110 R. 
Blackwell et al. [206] reported an observed value of 81 R at a rotational tem- 
perature of 294 K. The following year Krassovsky et al. [204] , from observa- 
tions at Yakutsk, USSR, reported a value of 140 R with a dependence on the 
rotational temperature. They also observed a seasonal variation with a winter 
maximum at Yakutsk. In the same paper they reported that observations at 
Zvenigorod, USSR, showed neither a temperature dependence nor a seasonal 
variation but a definite enhancement on cloudy nights. G. Kvifte [207] reported 
an intensity of 160 R for the OH (9-3) band. Huruhata [205] found that large 
variations in the intensity of hydroxyl emissions in one night were due to 
ascending and descending motions of the emitting layer. He also found that the 
intensity increased when the 10 mb pressure surface (approximately 30 km) 
decreased. However, a review paper by Krasovsky and Shefov [208] attributed 
the variations in intensity and rotational temperature to spatial patchiness and 
not vertical displacements. Greenspan [209] in an analysis of the spatial 
variability of the OH emissions reported that the latitudinal distribution showed 
a minimum near 30°N. G. J. Kvifte [210] reported a value of 90 to 110 R for 
the OH(9-3) band over an observatory in Norway, while Chiplonkar and Tillu 
[211] reported an intensity value of 113 R over Poona, India. L. L. Smith [212] 
reported that observations over Hawaii indicated that OH (9-3) band intensity did 
not covary with the 6300 k OI emission intensity, while Shefov [213] showed 
that there was a correlation between the OH intensity and the solar cycle. He 
also found that seasonal variations in wind speed and direction did affect the 
seasonal variations in OH intensity. Sivjee [144] found no enhancement of OH 
intensity during either electron or proton auroras. Thus, it appears that there 
are many inconsistencies and unknowns in the strength and variability of the OH 
emissions. From the sparse data available, Huntsville, Alabama, is apparently 
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located where there is a minimum in the latitudinal distribution of the intensity; 
however, this could be offset by the fact that there is a winter maximum in 
intensity of the upper vibrational levels [ 213] and also by the fact that the level 
of solar activity was relatively high on January 6, 1972. Observational data 
[144, 213] indicate a minimum in OH intensity between approximately 2200 and 
2345 CST. 

Ran and Kulkami [214] reported that contamination of the 6300 A OI 
emission intensity observations by nearby hydroxyl lines is very small where 
filters with half-widths < 10 A are used; however, the contamination could not 
be neglected for filters with half-widths > 40 A . They reported the following 
relative intensities for the lines: 


Line 

Band 

Relative Intensity 

6300 A 


3 

6264 A 


1 

6287. 6 A 

Pi(2) (9-3) 

0.0763 

6298.7 A 

P 2 (3) (9-3) 

0.0390 

6306. 8 A 

Pi(3) (9-3) 

0.0922 


This table shows that for every three photons of 6300 A OI emission there is 
one 6364 A photon, 0.0763 Pj(2) photons, 0.0390 P2(3) photons, and 0.0922 
P t (3) photons. 

The photometer used to obtain the observational data used in this analysis 
had a 6300 A filter with a 15 A half-width and an 80 percent transmission at the 
peak. Portions of the emissions in all three lines listed were transmitted 
through the filter and recorded as 6300 A OI intensity. 

The filter transmission function is combined with the above information 
on the position and relative intensities of these lines to compute the percentage 
of the observed intensity that should be attributed to the OH ( 9-3) band using 
the following: 
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I 


Total Relative Int. 


T 6287.6 R ^6287.6 + T 6306. 8^^6306. 8 

+ T 6298.7 R ^^6298.7 

(0.13) (0.0763) + (0.44) (0.0922) 

+ (0.73) (0.0390) 

(0.0099) + (0.0406) + (0.0285) = 0.0790 . 


The percent of the observed intensity due to the OH (9-3) band emission is, 
therefore, equal to the total relative intensity of the three OH lines divided by 
the relative intensity of the 6300 A OI emission, or 


0.0790 

3.0 


x 100 


2. 63% 


The Fabry-Perot interferometer used in this analysis had a filter with a 
5. 2 A half-width and a 55 percent transmission at the peak. A recent report by 
Hernandez 5 contains the new determinations of the positions of the Pi (3) and 
P 2 (3) lines as 63 06 . 83 7 ± 0. 009 A and 6297. 99 2 ± 0. 004 A , respectively. 

Figure D-4 shows the 6300 A filter transmission function. The positions of the 
three nearby OH lines are indicated on the figure, which shows that only approxi- 
mately 3 percent of the Pi (3) line is transmitted through the filter, approximately 
36 percent of the P 2 (3) line, and less than 3 percent of Ihe Pi(2) line. 

With these new wavelengths and an etalon spacing of 1 cm, which results 
in a free spectral range of 0. 198 A for this Fabry-Perot interferometer, it is 
possible to calculate where the above OH emission lines lie with respect to the 
OI line. These calculations show that: 

1. The P 2 (3) line lies between 0.054 and 0.062 A from the OI line with 
the most probable position being 0. 060 A from the OI line. 

2. The Pj ( 3) line lies between 0. 003 and 0. 021 A from the OI line with 
the most probable position being 0. 007 A from the OI line. 

3. The Pi(3) line lies 0.036 A from the OI line. 


5. Hernandez, op. cit. 
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WAVELENGTH (A) 


Figure D-4. 6300 A filter transmission function. 


Figure D-5 shows the position of these OH lines on a plot of an OI line 
shape which is a convolution of 18 emission line profile measurements on the 
night of January 6, 1972, between 1800 and 2330 CST. The average tempera- 
ture derived from this convolution is 841 K. From the plot it is impossible to 
ascertain if there is any emission from the Pj ( 3) line mixed with the OI line; 
however, neither the P 2 (3) nor the Pi (2) lines appear to be present. During 
this time period the OI line intensity decreased from approximately 100 R to 
approximately 18 R; therefore, it is believed that a P 2 (3) line intensity of a 
constant 2 R would be discernible on this plot. If this value can be used as an 
upper limit on the intensity of the P 2 (3) line, the results of Rao and Kulkarni 
[214] and those of Hernandez indicate that the intensity of the Pi (3) line should 
have been approximately 4. 6 R. The transmission function of the filter is such 
that only 3 percent of the intensity of this line is transmitted; therefore, only 
approximately 0. 138 R could be attributed to the OH emission. 

The individual emission line shape measurements for this night are 
shown in Figures 6 through 17. There is no indication in any of these measure- 
ments that the P 2 (3) emission line is present; therefore, it is believed that the 
contamination of the OI emission line by the Pi (3) OH emission line is 
negligible on this night. 
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ETALON PRESSURE 

Figure D-5. Convolution of 18 emission line profiles from January 6, 1972, 

with expected positions of OH lines. 


To support this contention, the data were subjected to two additional 
analyses. First, the data were input to a computerized version of a statistical 
analysis technique capable of separating a mixture of distribution into its 
components. This analysis showed that the data were not a mixture of signals 
but a single distribution. 

In the second analysis, several theoretical emission line profiles were 
generated and reduced with the results shown on Figure D-6. The values of the 
coefficients of the transforms for the theoretical 200 K profile and the 1000 K 
profile are shown as squares and circles, respectively. For comparison, the 
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NORMALIZED COEFFICIENT VALUES 





triangles show values of the coefficients obtained from theoretically calculated 
emission line profiles containing varying ratios of the intensities from a 1000 K 
and a 200 K source. The transforms from the observational data are also plotted 
on this figure with all values normalized to the value of the coefficient of the 
first transform of the 1000 K profile. Neither the values of the coefficients 
obtained from the analysis of the theoretical mixtures of sources at two different 
temperatures nor their trends with transform number have characteristics that 
are significantly different from those of a single temperature source. There- 
fore, this technique cannot be used to distinguish emission lines due to a mixture 
at two different temperatures from emission lines due to a single temperature 
source. 


The decrease in the observed intensity at conjugate point sundown pro- 
vides additional evidence that the major portion of the observed emission is due 
to the atomic oxygen transition, since the OH emission cannot be affected in any 
way by conjugate point conditions. 

For these reasons it is believed that the observational data used in this 
analysis are from the 6300 A OI line uncontaminated by any 9-3 OH band 
emission. 
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